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FOREWORD 


The domain of statistical physics has extended from the traditional 
problems in quantum mechanics, the solid states and the kinetic theory of 
gases to other nonlinear prdbl&ns in hydrodynamics, plasmas and optical 
syst&ns, for which we need to consider the fluctuation effects. Concep- 
tually, it has been found very helpful to view these statistical probl&ns 
frcrni a basic physical point of view which emphasizes their structural 
similarity. 

The study of stochastic systems has developed in two directions. First, 
the transition from a laminar state to a chaotic or turbulent state through a 
sequence of bifurcations shows a certain universality (see "Chaos and Univer- 
sality, " Nordita Selection, 1981). Secondly, when the system has reached a 
state of fully developed chaos or strong turbulence, the statistical methods 
used for their theoretical analysis find again a certain universality. A 
proper understanding of the many basic phencmena in astrophysics, space, 
atmospherical and optical applications depends critically on our ability to 
analyze the turbulent characteristics and the collective processes in these 
nonlinear systems . The lack of a suitable methods of treatment and the dif- 
ficulties encountered even for the simplest form of incompressible, homo- 
geneous and isotropic turbulence, have hindered the theoretical development 
of strong turbulence. 

From the physical point of view, a turbulent state is characterized by 
its transport properties: eddy diffusivity, eddy viscosity, coefficient of 

damping, or amplification. Their analytical determination requires a trans- 
port theory. The kinetic method is best suited for this purpose. The edcfy 
transport coefficients, as induced by the fluctuations of small scales, govern 
the evolution of larger scales. This requir&nent of scaling leads to the con- 
cept of renormalization groups, from which we develop the group-kinetic method 
of turbulence. We have applied the method to problems of atmospheric turbulence. 

The group-kinetic method combines the advantages of the kinetic method and 
the group-scaling. The kinetic description has the advantage of transforming 
the system of hydrodynamical equations that govern a nonlinear stochastic system 
into a master equation of lesser nonlinearity. The group— scaling enables the 
determination of the transport properties and the spectral structure by the 
one— point distribution function alone, without involving the two-point distri- 
bution. Our closure is obtained by a memory loss in the relaxation process and 
not by truncation of the infinite hierarchy of n— point distribution functions or 
correlation functions. 

The six verbatim sections include research accomplished during the second 
year, September 1, 1982 through August 31, 1983, of this two-year contractual 
effort . 
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SECTION 1 


General Considerations on the Group— Kinetic Theory of Turbulence 

I- PRINCIPLE OF THE GROUP-KINETIC THEORY 

The Navier-Stokes equation of motion with externai and seif- consistent 
forces is transformed into a master equation 

A. Z' 

+ L)f(t,x,v> = 0, 

/N 

for the distribution function f(t/X,v) of veiocity v ,with the equivaience 
reiation 

f(t,x,v> =P o jy ~ u(t^x)J 3 

that is vaiid for a constant or variabie density f ^ and for a fiuctuating 

/S /> 

fiuid veiocity u(t,x). Here L is the differentiai operator for the perturbation 
of the trajectory in the phase space. 

The totai distribution function 

f = f + f 

- / V 

consists of a mean distribution f = \ f / and a fiuctuation f . The usuai 
Fourier decomposition of a fiuctuating function contains too many minute detaiis, 
and a coarse- graining procedure is necessary. In analogy with the "renormalization 
groups", we decompose f into groups 

f° , f 3 f" 

of decreasing coherence, representative of the spectral evolution, eddy diffusivity, 
and relaxation, respectively. By formulating the relaxation as a functional of 
the diffusivity, we obtain a closure. 

The group- scaling has the advantage of determining the spectrum from the 
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singlet distribution function alone, without involving the 2-point distribution 
function as is required in the conventional methods of statistical mechanics. 

The kinetic equation is derived for . it takes the generalized form 
of the 'Fokker-Planck equation, with the diffusivity j- as an integral 
operator to represent the memory and the collective behavior. It is converted 
into the hydrodynamic equation in the group form, from which we determine the 
eddy viscosity 

K* r F(k)/ parameters 3 

vS 

and calculate the spectral function F(k). The parameters include the Coriolis 
field and the length of stability from shear and buoyancy. 

II. SPECTRAL DISTRIBUTION OF TURBULENCE 

A The application of the group- kinetic method to the Navier-Stokes 
equation of motion finds the direct and reverse cascades for the transfer across 

a spectrum, and derives the spectral distributions in inertia (k shear 

“1 —3 — 4 . 

<k and geostrophic turbulence (k and k 0. In boundary layer turbulence, 

the spectra of velocity fluctuations in separate directions can be treated, 
and a spectral gap in the stable boundary layer is expected. 

B. The extensions to the Zakharov equations, the nonlinear Schrodinger 
equation and the Korteweg-deVries equation are important for analyzing the new 
dynamical characteristics and transport properties of soliton turbulence in fluids, 
plasmas and solid, and the self- focusing in optical turbulence. 

III. NUMERICAL MODELING 

The kinetic equation for f and its conversion into fluid representation 
forms the basis for the modeling and prediction of turbulent profiles and 
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spectra in the turbulent boundary layer. Here the transport coefficients become 
explicit functions of t, x and the governing parameters, once the spectral 
laws are analytically determined by our group- kinetic theory. The pressure- 
velocity correlation Is also analytically derived. 

Our first aim in the modeling of turbulence is to provide an analytical 
basis for the Monin-Obukhoff similarity by determining the structure of the 
universal functions. Then, v\e can devise a pew numerical modeling which will 
incorporates the effects of scales, the collective behavior and the memory. 

IV. ATMOSPHERIC BOUNDARY LAYER 

Vi can divide the turbulent atmospheric boundary layer into: 

A. The lower boundary layer, or surface layer 

B The upper boundary layer which includes the free- convection layer, 
the mixed layer, and the Ekman layer. 

surface layer extends from the ground up to a height of less than 100 meters, 
and is characterized by the constant fluxes of momentum, temperature and 
humidity. The main tasks are the prediction of profiles and spectra by analytically 
determining the universal functions that could not be determined by the errpirical 
similarity theory of Monin and Obukhoff. 

The upper layer contains the additional effect of the Coriolis force, 
which enters into the differential operator in perturbing the trajectory. 

The above division of the atmospheric turbulence into two parts, as 
dependent on the absence o.r presence of the Coriolis field, is analogous to 
the classification of plasma turbulence into two parts: the Langmuir 
turbulence (i.e. without magnetic field) and the magnetized plasma turbulence. 
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SECTION 2 


Grotp-Kinstlc Theory of Two-Dimensional Oeostrc^fiic Turbulence 

C. M. Tchen 

The City College and The Graduate Center 
of The aty University of New York, N. U. 10031 


ABSTRACT 

The two-dimensional geostrc^fiic tuzbulenae driven by a random 

force is investigated. On the basis of the Idouville equaticxi ibich 

simulates the primitive hydrodynamical equations, we develop a group- 

kinetic theory of tuzbulence and derive the kinetic equation of the 

scaled singlet distribution. This distribution will suffice for the 

investigation of the spectnm of tuibulence, without having to resort 

to the pair-distribntion as tras with the usual kinetic theories. The 

collision integral has a memory and describes the pair interaction and 

its enhancement by the multiple interaction. 

Our kinetic equation of turbulence is transformed into an 

equation of spectral balance in the equilibriim and non-equilibrium 

states, as governed by the direct cascade and the reverse cascade, 

-3 -4 

respectively. The sequence of power laws k , k of velocity fluctuations 
are derived. 
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I. INTRCDUCTICN 

The two-dimensional geostrc^^c turtvlenae has attracted much attention 
[pj^rtoft, 1953; Chancy, 1971; Rhines, 1979] . Theories based ipon the 
diffusion approximation [ Leith, 1968] , the spectral properties llbvikcv 

9 

1979 ], the direct interaction approximation [ Kraichnan, 1967, 1971a, b ; 
Bouquet et al., 1975; Salmon et al., 197S] , and numerical modeling and 

simulations f Lilly, 1969a, b, Basdevant et al., 197S] have found a spectral 

-3 -4 

k ■ This spectrum changes into k if the geostrophic tuibulenoe is 

driven by a random force [ Saffman, 1971; Thompson, 1973 j . Deviations from 

these laws have also been discussed [ Gage, 1979] . The sequence of appearance 
-3 -4 

of the spectra k , k , and the dynamics of the direct and reverse 
casca d e were analytically vague lack of a generalized statistical 

theory of the geostrophic turbnjlence driven by a random force. 

Our investigation is divided into two parts. The first part develops a 
grovp-kinetic theory on the basis of our earlier scaling procedure £ Tdhen, 1978, 
1979 ] , and derives the kinetic equation of turbulence. Usually the analysis 
of the specbral structure requires either a pair-distribution function or a 
detailed Fourier decomposition. These are not needed here, because of the 
scaling procedure mentioned. Our kinetic approach has the advantage of not 
only including the random force in our homogeneous Liouville equation, but also 
of describing the interaction between the wave and the fluid pajdzicle, as 
characteristical of the large scale turbulence. The collective collision 
obtained will include the memory, the pair-collision and the multiple 
collision [Tden and Mis gulch, 1982 J . Hie latter effect is essential and 
will be applied to the treatment of the reverse cascade that is often found 
with the large scale turbulence. 

With the kinetic foundation described above, the kinetic equation of 
turbulence is transformed into its hydrodynamical form, and is used to 
develop the specbral theory of geostrophic turbulence in the second part o f 
this paper. The two parts are separately self-consistent and can be read 
indep&id&itly. 
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11. MICRODYNAMICAL STATE OF TURBULENCE 


The Navier- Stokes equation of motion 


(■^ ^ + u. 7 7^) « 


E 


( 1 ) 


of an incompressible fluid satisfies the condition 


V' u = 0 . 


( 2 ) 


% taking the curl, we obtain the vorticity 




(3) 


and transform (l)into the following vorticity equation of the two-dimensional 
geostrophic turbulencei 


+ u-7 - 

t ^ ^ 


(4) 


Here we have 


'n “ (Uj^, U2» 0) and ^^=(0, 0,^). 


(5) 


The field 


^ a. A 

E . Ep E^ 


( 6 ) 


has two components. The component 

EpS-^7P . 

J ^ 

which is the gradient of pressure “p at a constant density ^ j is not 
present in the vorticity equation (3) , because 

H7xEp ®--7xtTp = 0 . 


(7) 


C8) 
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/\ >N 

The other component E forms a finite vorticity source '^X E . 

The equations (1) , (2) and (4) describe the microdynamical state of 
turbulence, and are used as the point of departure of our statistical 
treatment. In our transport theory by eddies larger than the viscous cutoff, 
the effect of the kinematic viscosity will be neglected, but will be 
restituted in the spectral balance where the viscosity gives a molecular 
dissipation. 

For the development of a kinetic method, we write the microkinetic 
equation in the form of the Liouville equation: 

+ L(t)J = 0, (9) 

with the differential operator 

L(t) » V V + ' (10) 

The detailed distribution function f is normalized to unity, as 

fdv f(t,x,v) = 1, (11) 

and assumes the form 


f(t,x,v) 


^^-^(t,x)J , 


( 12 ) 


in order to be consistent with the Navier-Stokes equation (1) and the 
equation of continuity (2). It is not difficult to show by taking the first 
two moments j that Liouville equation will reproduce the two hydrodynamical 

equations Cl) and (2). The Liouville equation is homogeneous and has less 
nonlinear terms than does the Navier-Stokes equation, because the function 
u(t,x) of velocity of fluid is replaced by an independent variable v . 

We retain the only nonlinear term E '^f in order to describe the 
mode- couplings. 
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III. GROUP-SCALING PROCEDURE 


The Liouville equation (9) contains too many minute details which are 
unncessary for a statistical treatment of turhulence. We apply a scaling 
procedure hy dividing the detailed distribution 

? = f + f (13) 

— <0^ 

into a mean distribution f and a fluctuation f , and the fluctuation 

f « f° + f (14) 


into a macro-group and a micro-group f ■ The function f(t,x,v) 

gives the distribution of velocity in non- equilibrium. The macro-group 
calculates the spectrum of turbulence'without the need of developing a 
separate kinetic equation of the pair- distribution function <? . In the 

following, we develop a kinetic equation of the macro-group and 

investigate the spectrum of turbulence. 

The evolution of the macro-group is controlled by the transport 
property as shaped by the fluctuations of the micro-group, while the 
approach of the transport property to equilibrium is obtained by the 
relaxation. These three transport processes of evolution, transport 
property and relaxation are represented by the three groups 

f , f" 


of the distribution function, or by the three groups 

e", E' , E" 


of the field. The interactions between the groups are described by a 
coupled system of equations. 

The groups have their time scales characterized by the correlation 


times 




T? 

C 


It 


(15a) 
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For formulating the transport coefficients, the correlations, as given by 
the groups E' (t-ZT ), E”(X-X^), are integrated with respect to 'V and 

T; ' in the domains (0, t) , (0,X ), respectively, so that the evolution 
times t ,T>T 'can be ranged among the correlation times (15a) , in the 
following manner : 


t?T/ > x'' T (15b) 

The degradation of coherence (15) will constitute a property of 
quasi-stationarity of one group with respect to the other. 

The individual fluctuating groups 


can be decoirposed into Fourier components with overlapping wavenumbers . But the 
global averages > e g ^ 


are deterministic, and are separated by their adjacent wavenumber domains 

(16) 

(17a) 
(17b) 


(O.k) , (k«fl) , (k'^k.bo) 

with k, k* varying from 0 to^O . 

For the sake of convenience, we introduce the scaling operators 

A , A° , A' , A" 

- A + A°, 1 - A, A^^^=A' - A" . 
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IV. SCALING OF THE LIOWILLE EQUATIONS 


By applying the operators A and A' to scale the Liouville equation 
(9), we obtain the system of coupled equations: 


ca,. + A°L)f° = - L°f + C° , with C°s-A°L'f 


(■b ^+L)f 


- L'f^ - C . 
o o 


The second equation of the system may he equivalently replaced hy 
0>^+A'L)f = - L'f^ . 

Upon integrating (19) and (20) , we have 


"A* / dXUCt.t 

^0 

■ *7: 


-l) L'(t-T) f^(t-T) 


and 


f » 


dr U(t,t-T) C^(t-T), 




L 


( 18 ) 

( 19 ) 


( 20 ) 


( 21 ) 


( 22 ) 


.. '0 

Here U and/l, are evolution operators as related to the differential 
operators L and A'L, respectively. The initial value f (0) has been 
neglected, since it cannot produce a finite correlation A^L'(t)f(0) 
at large t , by (15). 

Subsequently, we prernultiply (21) and (22) by -A°L' to find the 
collision in the following two forms: 


c” - (p 




dT L'(t) U'(t,t-T) C^(t-X) 


C° " A° C , with H° s -A°L' (t)U' (t,t-T) 


( 23 ) 
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C" =A° 




(24) 


The collective collision G is seen to consist of tro parts: The pair 
collision 


(C®= k° 


( 25 ) 


from the correlation between two micro-fields a'^L ' A* UL ' , is proportional 
to a' • The multiple collision, from the correlation A*^L ' U ' between 

a micro-field and the cluster of the organized micro- fluctuations, is 
proportional to ^ ■ 

The governing collision operators 

rt 


A°A' = A° / dC L' (t) A'U(t, (i-zy 

0 

. ■i.A’P'-'i 

A°A'~ A° fdx L'(t)A'A(t,t-X)L'(t-T) 


(26a) 


(26b) 


are defined by the diffusivities 


A°jp' «= A° fdT K'(t)A'U(t,t-t. )&’ (t-^ ) 

A'^D' = A° dV E' (t)A'A (t,t-T)12'(t-C ) 

«H / •«» 

''o 

which themselves m^ serve as operators. In (23) we have replaced 


(27a) 

(27b) 


A°L*(t)A'U(t,t-T)C^(t-x: ) by A°L* (t)U* (t,t-i;)C^(t-r ) 


without loss of g&ierality. 





2-9 


The operator U(t,t-i:), with t' » t-T , is governed hy the Liouville 
equations: 

+L(t)] U(t.t’) - 0 (28a) 

- L(t')j U(t,t') = 0 . (28b) 

The operator A is related to U , and the relation has been developed 
by Weinstock [l969j . 

By scaling (28a) by means of A" , we obtain the kinetic equation 

( + L^)U° = - L°iT + H° , (29) 

with the collisions 



H°(t,t’) = A° A'|u°(t,t»)j 

(30a) 


If (t,t«) - , 

(30b) 

as related by 





(30c) 

i n analogy with 

(23) - (25). 
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V. PPOPAGRTOR AND PPOBRBILITY OF RETROCSiADE TRANSITIOf 


-A 

ZTie Liouville equations (9) and (28) which govern f and U have 
the same differential operator ^ , and therefore the same characteristic 
equations governing the dynamical variables in the form 


dx(t') ^ dC(t') 

= v(t') , — =• E(t’) , (31a) 

dt' '■ dt» 

or equivalently, 

d^x(t’) 

^ ■= E(t') , (31b) 

dt'*^ 


with the conditions 

x(t) * X, v(t) « V, (31c) 

by the nature of the differential equation of the second degree, ndiere 

X and V are two independent variables. The values (3 Ic) will be called 
Idle "initial conditions " in the retrograde transition from t to t ' . 

By a change of variables 

*S. 

xCt-"^) = x~ y^(r) , y(C-t) = V “ V(t) f- , (32a) 

^ MW 


we can write the characteristic equations in the integral form 

XiT) =fdrv(t-r) 

Jq "" 

v(-c ) “f dr E(t-T ) . 

■0 

Note that x(t-r) from (32a) can be rewritten as 

x(t-T ) >= x"/* t-T^ ,K~Jl (-0 ( t,x)~l 


(32b) 

(32c) 
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or as 


x(t--C) 



in the Lagrangian representation, with 

£ ('D 1 t,x) replaced by J 

for the sake of brevity. The vertical har in I denotes that the 
Lagrangian displacement J^(x) was made during a time interval X along 
the trajectory that passes hy the point x at the time instant t. 

The systems (31) and (32) are called the equations of the path 
dynamics in the differential form and the integral form, respectively. 
Both forms equally describe the trajectory of a fluid particle which 
occupies the position x(t') at the time instant t' , provided 
the trajectory passes by the point x at the instant t while having 
a velocity v . 

It is to be stipulated that the essential function of the propagator 
U IS to impose a Lagrangian representation of the function, say E as 
in (27a), in the form 


Eft',K(t’)] S U(t,t’)l(t'), (33a) 

following the trajectory that is determined by the differential path 
dynamics (31) . 

Alternatively, with the path dynamics (32) , the Lagrangian 
representation can be obtained in the integral form 

E |Tt-T J ' di E(t-T,x-J^ ) '^(t-T,x-^jt,x), (33b) 


by using the probability density 

'p(t-r, x-£ j t,x) . (34) 

«* 

The integral relation is written for the purpose of transforming the 
Lagrangian function into the Eulerian function. 
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The probability density satisfies the condition of normalization 


The density (34a) 



1 . 


defines the probability 


p(t-r ,x-i 


t,x) di, 


for the fluid particle to occupies a position between 


x-^ and X -(i +dH ) 


(35) 


at the time instant t-T , provided it follows the trajectory which 
passes by the point x at the time t. For the retrograde transition to be 

specified by the integral path dynamics of ZMt as given by (32), we 

T '*■ 

write [Tchen, 1944j 

p(t-T, x-X 1 tjX) = ]^(t-T, I t,x)l. 


or briefly ^ 

''pC-r.Jl) - • 

^ w* -V ** 

The abbreviated form indicates a quasi- stationarity of p in the t,x 
space, so that the Liouville equation can be written in the form: 


+^(t-r )l 'p(t-^,x-4 I t,x) = 0, (37a) 

si 

with the differential operator 

Xf(t-T) = -<^(t-T).^ , 

It is not difficult to show by moments that (37a), together with (36), 

will reproduce the integral path dynamics with the initial velocity 

V that is a random variable having a distribution f(t,x,v) . By 
^ /\ *** **** 
definition (12), v can be identified as u(t,x) in the return to the 
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micro- hydrodynamical description. Thus the Lagrangian representation (33a) 
hy means of the propagator U(t,t') is straight forwardly kinetic. Qi the 
other hand, the Lagrangian representation (33b) hy means of the prohahility 
P(T.P) i s valid in both the fluid and kinetic descriptions. In the former, 
the initial velocity is u(t,x), and the Liouville equation (37a) becomes 
dissociated from the distribution f(t,x,v). In the latter, although the 
initial velocity v is not directly involved as a variable in the 
probability density ^p(f ,ji ) , it is implicit in the path dynamics through 

the distribution f(t,x,v), so that the same Liouville equation can be 
copsidered as an equation of a kinetic significance in the contracted 
dimensionality, i,e. in the sense that the governing path dynamics requires 
an initial velocity to be determined kinetically. The contracted dimensionality 
gives to the Liouville equation a simpler differential operator that is 
independent of Jl and therefore a simpler Fourier transformation. 


In the following, we shall exploit this advantage of the contracted 
dimensionality by using the probability of transition in transforming a 
Lagrangian function, e.g. (33a), into an Eulerian function, e.g. (33b). 


J 



2-14 


VI. MEMDRY m THE COLLECTIVE COLLISICN 

The problem of the derivation of the collective collision C° in 
terms of the pair collision involves the solution of the covpled 
integral equations ^ and (30) , and the treatment of the Lagrangian 
functions associated with the propagator and its collision, from (29) 
and (30b) . The problem is ccaplicated. [i/e shall make the simplification 
by separating the chain of correlations (15b) into two sections 

t 

and 

T 7 T/ > T ; (38b) 

referring to the processes of transport coefficient, i , e , diffusivity, 
and relaxation, i.e , propagator, respectively. The first section which 
has a correlation of long duration shall preserve its memory, vhile the 
correlation of the second section vhich has a short duration should not 
propagates its memory beyond the time spaa 'O and exert an effect on the 
evolution of f°(t*x,v) , 

With this memory-loss we simplify the equations o f collisions (2S) 
and (30b) into the system of equations 

C° -HI ^ C° 

H IF , 

which we can ccnbine into 

C° =C ^ C° . (40) 

Jt leaves the Lagrangian-Eulerian transformation by the method of the 
transition probability, as follows: 

C° U * C° 

Jdvfdp(r,h C°(t-T,x-i). (41) 

✓ o 


(39a) 

(39b) 


i 



W see the transformation 


U(t,t-D C°(t-T) “fdl p(t,/) C°(t-T,x-i) 

from (33b), in which the kinetics of U as governed hy (29) is replaced 
by the dynamics of "p" as governed hy (37a) . Note that 

C° H , u* 

are in the t,x,v space. Qi the other hand, in the "T , space, P(V,H ) 

«« -A ^ 

does not carry v , except in the path dynamics as an initial condition. 

•V* 

Further simplification is obtained hy a Fourier transformation of 
(41) into 

(27t)"^ p(w,k) C C°(^,k), (42a) 


with 

0^ («).k) 

in d-dimensions, 

form 


CO 


“ -^'>Ygtei(<J-3s:.«)re-ec,k) (42b) 

so that t^^e integral equation (40) is transformed into the 


c°(^d,k) = C\cJ,k) +oe 0,k) c°(u) ,k) , 


or, hy collecting C*^, 

C°(cJ,k) = . 

W / 

with 

^(co.k) = [ 1 -0^(<J,k) 1"^ . 

The prohahility function' PC"C,| ) is determined hy the following 
differential equation of evolution: 




I 


(43) 

(44) 


(45) 


that is scaled from the Liouville equation 
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with a fluctuating differential operator. We find the solution: 

r*' 

P(T.k) = exp [- , 

0 

Frctn the path dynamics (32) , we find 

= T ^ [ <D)°(r)) + (c)} ] , 


(47) 


(48) 


and that(^' (r)^ is asymptotic and ^1D°(X)^ is not asymptotic in the 
evolution of P(T»1) from (45) , 

T - 

Now we can calculate!^ from (42b) , ty substituting for P fran (47). 
Wfe find tkat 


with 


[Si (y,k) = c W ff(£4,k) , 

o< 3> 


Di C< -f ^ 


is controlled by the memory— loss function 

tCO 




V A 


^ la-p -4,r"-V^r 
qr T -e -e- 


which determines the life-time of memory in the form 


T ^ 
H 


'/i 


The characteristic frequencies are 


k'V , 0^ , m° (c° k^) 




and the numerical coefficients are 


C =7T/c’ 
Of 


c=(l/3)^^^ 




(49a) 


(49b) 


(50a) 


(50b) 


5Ial 


(51b) 
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As seen from (43) and (44), the collective collision C is 
related to the pair collision hy a factory^ , which can he revritten 


as 


ri 


= [(/- «,/ ^ <J , 


(51) 


in terms of the real and imaginary parts of (31^ , hy (49). Hence the 

determination of (51) rests on the calculation of the memory integral 
(50a) and its decomposition into its real and imaginary parts. The 
calculation is made hy an interpolation over the three regions dominated 
by 4 » and m° , separately. The results are given as follows: 

(a) For weak turbulence, i.e. , we find a shielding 


in the collision hy a factor 

/ 2 = /- - I ■ 


(52) 


(h) For strong turbulence, l.e. fl«. oV, , we find an 


enhancement in the collision by a factor 


Jt> 


if m ° << 


^ 12.51, 

[l + 2.77 (4>'/m°)^^ 1, if . 


'■p-- ’ y 

W can approximate ^ by its constant asymptotic values which are ; 

(a) weak turbulence , Cl» / In® 

P 


(53a) 

(53b) 


(3 ^ . 

(b) strong turbulence. 




(large k) 


^ ^ ri, for 

^ 2.51, for ^ ;^^ij*(sinall k) 


(54a) 

(54b) 

(54c) 


3 
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Xt entails the siirplification 

C°(t,x,v) ^ (C (t.x.v) 

of (43') . 


(55) 


VII. SPECTRAL BALANCE 

The kinetic eqation 08) of , rewritten as 

= - L°f + C'^ (56) 

has a collective collision vhich can be written as 

c° . (57) 

with a diffasivity found by (55b) , 

Consider now a homogeneous tuxhulence, i , e , with ~u 0, E 0. 

By taking the first moment of the kinetic equation ( 56 ) , we obtain the 
following equation of macro-ncmentxMi: 

i'h + A° u° P7^)u° - E° + Q° , (58) 

with a collision 

by (57) . Subsequently , -we take the curl of (58) and (59) , to obtain the 
following equation of macro-vorticity of the two-dimensional geostroffiic 
tuzbulerxoe: 

( > . + u° 7 -V 7^)t° - V k E° ^ Vx (60) 

^ -V -* F vm ^ 

<w% 

Finally, ipon multiplying (60) by and averaging, we find the equation of 
spectral balance in the form: ”*** 
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( 61 ) 


with the following transport functions: 


cascade tranfer 


t : = <■?" 

molecular dissipation £ ~ 


(62a) 

(62b) 


Note that the field 


E° - E ° 

^ «<* A. 


consists of two parts: The first part has a zero vorticity hy (8), and 
the second part, which is a random source, does not correlate with 
Thus the field E° does not appear in the spectral balance in the form of 
the correlation . But it governs the collision 

Q°= Qp° + Qx° * 


with 




through the diffusivities 

<3)'> -<d;> 


(63b) 


(63c) 


with^ ^ -t£) 

<]pp^“ J dt/Ep' (t,x>4t^,t-0£p’ Jdtr<4^' (t,x)y^U^,t-mx’ (t-Tjj> .(63d) 

0 ^ ^ 0 ^ 

By separating the collision into two parts as in (63a) ^ we find the 

transfer function also separated into two parts, as follows: 


+ T. 


X 5 


with 


Tpo 7X Qp7 » V 7aQx° ) * 

In the following, we shall follow the same method of treatment of 


(64a) 

(64b) 


J 
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the transfer function as developed earlier by Tchen fl98lj . As n otations 
we introduce the spectral functions ^ (k) , and the spectral 

density , such that 

(n°^} - 2 rL' zAL'F (k'). if (65a) 

J n ^ V -w 


0 

and the second spectral moments 

k 


R ° B 2 f 

r 


,.2 


dk' k'*- F (k’), R. 

0 r 




dk' k'^Fy(k'), 


(65b) 


with the relationship 


F^(k) =■ d"^ k^F^(k), (65c) 

by definition (3) in d=2 dimensions. omitting the details of calculation, 
we find the results as follows: 

(a) Cascade-transfer by the E „ “ fluctuations 
The transfer function is obtained in the form: 




for large k 



for small 


k . 


(66a) 

(66b) 


The coefficient ^ represents the effect of the multiple collision (54). 
The transfer function is seen to be governed by the eddy viscosity 



which has an asymptotic value 



(F^(k")/k"^ ) 


1/2 


for large k , and by the rate of damping 


(67) 


J 
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A' 

with oD 

A \ / V r Hit" f (k“) 1 

<«<“'> J 

for small It ■ Note that (69a) is the solution of 

The asymptotic formulas (67) and (69a) agree with our earlier cascade 
theory [Tchen, 1973, 1978 J and other dimensional theories 
^Heisenberg, 1948; Gisina, 1969^ . 

(h) Cascade-transfer hy E„ “ fluctuations 

The transfer function is obtained in the form 


''' \ o 

jj/R , for large k 


, for small k , 




and is governed hy the eddy viscosity 


<sj^0^')) G(^') 


for large and the rate coefficient of damping 


ft " » [ dk" k"^ (k"^ G(k") 

^ \J ^ <V 

for small k . Here the modulation function 
is characterized hy the two frequencies as fqJ.lows! 
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It is to be noted that the simplified formulas of all the transfer 
functions listed above are characteristically in the form of products of 
two functions. A general formula not in the form of a product has also 
been developed by Tchen (1981b), but is too complicated for the application 
to the present problem. 

We first consider the portion of the spectrum of large k . The small 
eddies are embedded in a gradient of big eddies and in their coupling with 
the bigger eddies, the small ones play the role of an eddy viscosity, 
to cause a loss of the bigger ones at the rate . This transfer 

is in the pattern of a cascade from the big eddies to the smaller ones. 

It is of the gradient type and is called the direct cascade. New we 
consider the portion of the spectrum of small k ■ There the big eddies have 
no more gradient left in the medium in which they evolue, if the turbulent 
medium is isotropic. The coupling between the big and small eddies i s 
not through the gradient transfer, but is of the type of the damping. The 

rate coefficients n' - n.,' are originated from the wave-particle interaction 

.. ^ A 

[Tchen, 1981 J , and may become negative when k is sufficiently small. 

Then they give an amplification. What this happens, thb transfer functions 
reverse their roles from a loss into a production, i,e, > o . 

Such a transfer mechanism is called the reverse cascade ■ An analogous 
phenomenon is known in plasmas as the Landau damping or amplification. 

If the demarcation wavenumber between the direct and reverse cascades 
is ^j.gy 5 the large and small values of k refer to 

k > k and k k , (75) 

^ rev ' rev 

respectively. 


1 
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VIII. STATISTICAL EQUILIBRIUM AIIB NCN-EQJILIBRIUM 


By differentiating the equation of spectral balance (6L) with respect 
to k, we obtain: 


F^(k) = - - T ° , 6 . (76) 

*^5 K ^ 

The left hand side vanishes in a statistical equilibriim, yielding 


t ° + T° = £ . (77) 

The constant of integration i n the ri^t hand side is determined by the 
condition at k ^ oo , i . e , 




and 


In addition, we note that 


__e» 


(78) 


T^ ^ 0, at k=0 . 


( 79 ) 


For the inertial subrange in statistical equilibrivm, we can write 
the spectral balance in the form 


T T° + T° 

ip -r i^ 


(80) 


by emitting the viscous dissipation t ° by definition of the inertia 
subrange. Here is a sink in the enstrephy transfer across the spectrum. 

The statistical equilibrium requires that the two transfer functions 
find a net positive transfer to balance . 

On the other hand, if the net transfer is negative, ^ ceases to be 
a useful parameter, and the hypothesis of the statistical equilibrium 
becomes invalid. Consequently, we have to return to the original equation 
(6L) for the spectral balance in ncxi-equilibriim: 

- 1 ^ ( 81 ) 

by assuming a supply to the spectrum from an non-stationary source of larger scale, 


1 
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Here we have again omitted the molecular dissipation. The negative net 
transfer, i.e. amplification of the enstrophy in time, 

and forms a reverse cascade toward the small wavenumbers spectrum. 


IX. SPECTRAL DlSTRIBUnONS 


In the preceding Section, we have distinguished between a direct 
cascade and a reverse cascade, as characterized by a transfer of the gradient 
type at large k and by a transfer of the damping type at small k . W 

analyse the spectral distributions for these two cascades in the following 
lines. 


A. Direct Cascade 

W consider the joint enstrophy transfer by u and S “ fluctuations 

A 

in the spectral balance (80), and rewrite it as 

'-k. 




1 j dk’ k'^ F (k') = ^ , 

Jq ^ > 


(82) 


with the use of (66a), (67), (70a) and (71). The eddy viscosity ^ 
governs the transfer under the driving force exclusively, and the eddy 
viscosity governs the transfer without the driving force. 

An approximate solution of (82) can be obtained by an interpolation 
of the following two separate equations of spectral balance: 

A 


;(|K^^^Jdk* k'^ F^(k') 




'k j 

dk’ k' F (k’) 




(83a) 

(83b) 


finding the asymptotic solution 


F (k) = C k"^, C « 2.6 

It ' ' ^ ' 


u 


(84a) 


of (83a) , with 


J 
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bjj f^/3 ^_2 ^ 2/C » 0.77 , (84b) 

and the asymptotic solution 

F^(k) “ k~^ , => Cj, ^ = 0.60 (85a) 

of (83b) with 

^®Sc* ^ ^ ^ ' (S5b) 


Note that, for obtaining the solution (85) of the equation of 
balance (83b), the eddy viscosity (71) has been written as 





G(k’) 


» <^s dk" 2^ k" G(k'')> for finite k, 

where is the trace of the spectral density tensor, i,e, 

<8 >' trace<SjJ ) . 


( 86 ) 


(87) 


As a white noise it has the form 


^^(k")^ “ (s ) , independent of k", for finite k" 
*0, for k’’“C0 


( 88 ) 


A/so note that 

dk" = 2TTk" dk" . (89) 

W have calculated G(k") from (73) and (74) for substitution into (86), to 
find^C^’^as was written in (85b). The numerical coefficient has been 
determined from the equation: 



( 90 ) 


1 
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Finally by substituting for thus obtained, we have found (85a) 

as the solution of (83b) . 

New we return to the problem of the interpolation of the two 
asymptotic solutions (84a) and (85a) , by incorporating the two eddy 
viscosities (84b) and (85b) in the equation (82) of the joint balance, 
which becomes : 

A 

"’k % 2 J JO F ft') - 

The equation (91) shows that the dominant eddy viscosities are^ti the form of 

for k < k^, and for kj> k^ (92) 

Then we find the following formula of interpolation: 

F^(k) - C f ^(k/k^) , 

with . ^ j . 

■ if ^ J 

'■(l+k/ko)2 1+ k/k^ 

and 

- aA„/c) {. 

~ 0.46 <^s^-l/4 

It is easy to verify that the general solut-on (93) include the 
asymptotic solutions (84a) and (85a) for k ^ k^ and k ^ k^> respectively. 
It is seen that the spectral law k~^ precedes the spectral law k”^ 
in the sequence of increasing wavenumbers. 

B. Reverse Cascade 

For the reverse cascade which appears in the region of small wavenumbera-, 


( 93 ) 

( 94 ) 

(95) 


3 
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we shall use the formulas (66b), (68) and (69), and write the spectral 
balance as 


or 






(96a) 


hy a change of notations from (65b) and (65c) . Here C is the rate of 
increase of 



(96b) 


(97) 


in non-equilihrium, and is assumed to he independent of kin the global 
balance which may include the supply from the synoptic scales. 

W divide isth sides of (96b) by and differentiate, to obtain 


t-K) = RJ 


(98) 


or 


where 


Rj = r/3'U ( /l; f ylj 


(99) 

( 100 ) 


is seen to be a positive and dimensionless number in this subrange, 
since * is a positive function which increases with k monotonically, 

Ud 

by definition (68b). We can estimate to be of the order of unity. 

Finally we differentiate (99) with respect to k and find the 
spectrum 


F (k) 
u' ' 




<s> 


-1/4, 


/.(k) 


-4 


V 




, -1 
K 


0.24 


( 101 ) 


] 


where 
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t - - 4 = (a:/cV‘/^ 

I + ¥i 


si: /a' 

is estimated to be 

a- - 2 

if we approximate /L' - been made. 

Since the reverse cascade occurs at small k, f\re have taken the value 
(54c) for ^ . 

The solution (101) takes the following asyirptotic expressions 


F„(10 - r C* , for k C 


k-‘>. k>k^. 


C*^f%{rr/,i)^^l.o'g^ 4Vo.^</ 


We note that the numerical coefficients (104) in the reverse cascade 
are smaller than their corresponding values (84a) and (85a) in the 
direct cascade. 


J 
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X. CONCLUSION 

In the direct cascade, a critical wavenmber separates the spectrum 

into tm regions: The region at Jt > k shows a dcadnant role of the 

-4 ° X 

random force, and has a power law k . Tne region at kx has a 
negligible effect of the random force, and yields a power law k ^ , In the 
reverse cascade, the same power laws repeat at a new critical wavenmber 
k©* ( ^ k©) , with their numerical coefficients modi fied ty the multiple 
collision , and with diffejrent parameters, A bmp spears in the transition 
from the direct cascade into the inverse cascade ^Weinstock, 1978J ■ 


J 
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EQOVALENr METHODS 

FOR QUASILINEAR TURBULENT TRAJECTORIES 


J. H. Misguich and C. M. Tchen* 


ABSTRACT. 


The propagator formalism is summarized in a self-consistent 
way for the asymptotic quasi-linear equation. 

A comparison is performed between the propagators and the 
Green 's functions in the case of the non— asymptotic quasi-linear 
equation. This allows to prove the equivalence of both kind of 
approximations used to describe perturbed trajectories of plasma 
turbulence. 


*This work investigates the detailed dynamics of the perturbed 
trajectory in turbulence. The manuscript is prepared for publication. 
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Renormalization introduced in the microscopic 
theories of plasma and fluid turbulence essentially consists 
in taking into account the turbulent perturbation in particle 
trajectories . 

These turbulent modification (perturbed orbits) can 
of course be calculated by integrating the equations of motion 
in the fluctuating field. In a kinetic theory one uses however 
more compact and powerful tools which allow to directly describe 
the time— evolution of observable averages of all dynamical 
function by means of a unique entity: the propagator acting on 
phase variables x and v, or the Green function. 


The aim of this work consists to prove the equivalence 
between two approximations which have been developed indepen- 
dently in these two foemaiisms, This allows ya to bridge the 
gap between different methods used by various authors. 

Turbulent modifications of the trajectories are 
described in the lowest approximation as quadratic function of 
the fluctuating field, and by considering all other trajectories 
as unperturbed by the turbulent field. This approximation has 
also been used in the quasi— linear equation for the distribution 
function and is thus referred to as the quasi— linear appro Klmat ion 
■for the trajectories { although it uses free trajectories as 
basic ones) . 


The turbulent trajectories so— obtained are then-nsed 
as a basic ingredient in the renormalized theory of plasma or 
fluid turbulence in the next approximation which is referred, to 
as the Renormalized Quasi— Linear approximation ( RQl{ 1 » ,5,6, / 
This latter appears to be equivalent to the Direct Interaction 
Approximation (DIA) introduced in fluid turbulence /i ,3/ . Thus 
the renormalized propagator appearing in this approximation ( in 
the weak-coupling limit) actually describes lowest order turbulent 
trajectories(^ 1 , e , in the quasi— linear approximation ) and 


i 
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we will limit ourselves to this approximation for the particle 
trajectories . 


About the methods used to describe these trajectories, 
we have to remark that one has first integratedthe equations of 
motion /!/. WEINSTOCK /2/ has then introduced in this probleti niathesnatical 
operators ' acting on the phased 'variables x and v , the propagators J7/ 

the role of which consists to describe the time— evolution of 
associated dynamical variables x(t) and v(t) , and, consequently, 
(Liouville's theorem) to describe the time-evolution of distri- 
bution functions. In order to calculate the action of these 
propagators, WEINSTOCK had however to introduce particle 
trajectories which have been calculated by integrating the 
equations of motion. Some subtle effects , like that giving 
rise to an average displacement of the particles have been first 
missed by this method /3/. 

Other unexpected effects, like the non— vanishing 
correlation between particle velocity and position, have been 
found due to the first explicit calculation of the turbulent 
propagator /8/. In Ref. 9 and 10 we have discussed the explicit 
relation between trajectories and the propagator for forward 
and backward propagation in time. 

Other authors have based the kinetic theory on the 
Green functions /6/ . We will compare here these methods and pro- 
ve the equivalence between the approximations used. In spite of 
the global equivalence a few differences remain, namely concerning 
non— Markovian or memory effects. 


J 
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IX. KINETIC QUASI-LINEAR EQUATION 


ll^A. Me th od of derivation 


The renormalized propagator and Green 's function 
which will be derived are described by the quasi-linear 
equation which will be deduced here in a simple way. 

For electrostatic plasma turbulence, as well as 
for Navier— Stokes fluid turbulence, the starting point can be 
put in the form of the plasma Klimontovich equation 

* V.l * CC£W'\ 1 

. "* — J 

( in the plasma case a ratio q^ is included in the electric 
field E_ , with the charge and m the mass of the particle) . 

In the plasma case H(X/V, t) is the Klimontovich microscopic 
distribution 

ti 

Here x and v are the phase space variables, while x^(t) and v^(t) 
are the dynamical (time-dependent) variables described by the 
exact equations of motion. The electric field is the microscopic 
one given in term of N by the Poisson equation /ll, 12/ , 


The ensemble average introduced in statistical 
mechanics allows one to define the usual distribution function 
by 




and the Vlasov mean field for inhomogeneous systems by 




average 


Fluctuations are defined as the difference with the 

f ' = N _ <N> S-'< 

e ' = E _ <e> J s- 
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The Klimontovich equation can be written in the contact form 


N = L H 


■jT.6 


where l_ s, i. -i- U* 

with _ . . V 

L=<L> = -v.V-<6<2,t)>-;^y 
L' r-E'Cx.t)-- 

- *T>y 

Separation in average and fluctuating equations gives 

-- LF' + L'f + -feL'P 


T.S 

J.9a 


I 

where = L'^' - <>-'f '> ^ Cl - L' F' 

A... > 

describes the fluctuating part of the product. Let us remark 
that this last equation ( TT4 ) also writes 

-iitF’r LF'h L'F -<L'F’> ^ 


From the two equations {l.$) andCn-^) a general 
kinetic equation can be derived in the following way : equation 
(X.Sj is solved for in terms of f and the solution is 
substituted into (l.S) ■ A closed equation is then obtained 

for f , which still contains the initial fluctuation^' (^c) i 
this equation is called the general master equation. This 
equation, obtained by Weinstock from the Vlasov equation, is 
fully analoguous to the master equation obtained by Prigogine- 
F.4 si bo is in statistical mechanics, or to the one obtained by 
Mo r i or Zwanzig . 

In order to go to a kinetic equation, irreversibility 
has to be introduced in some way : it usually consists to take 
the limit of long times compared to microscopic elementary 
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times {Euler ian correlation, time of fluctuating fields) and to 
neglect the term Involving the Initial fluctuation. 

One of us has developed with BALESCU a theory which 
shows that this general kinetic, equation can be obtained In 
a rigorous way for one ("kinetic") projection of the distribution 
function, due to the existence of two Independent subdynamics 
/13,14/. This theory brings justification and necessary validity 
conditions for the simple derivation which Is made here. 

II . B . Pro pa gators . 


The basic tool for solving equations (Xd*/bjis the 
propagator. We will first present the simple and trivial exenple 
of the free propagator associated with unperturbed particle 
trajectories. In such a simple case, the evolution equation 
reduces to 

\F = UF 


where 



V.-7 


The solution immediatiy writes 


Here the exponential operator 

Ut , -bv.Y 

e = = e. 

which Is solution of the same equation 

Uo (b,o) - L„ 


allows us to describe the time— evolution of the distribution 
function as 




} 
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It is inportant to remark that this same operator also 
describes particle trajectories: 

•JS v.^ 

Xtb4®) = ^ +^Y 

V Lh-*-a) r Uo Ct,\: + 5) Y. - ' "Y - Y 

and this allowsoneto check the "Liouville theorem" in its simple 
form 

^ (x,w,t)=e = f [_S-vt ,oJ 

=. F ,Y(-k),o] 


This is the basic advantage of the propagators: the same 
mathematical tool allows one to describe the forward propagation 
of the distribution function in time: 

fCx,v,\:+^)- UoCbi-t,t) Y.Yit) 

and the inverse (backward) evolution of the dynamical variables ; 
X(t)x Uo + Xd-iZ.) 

thus 

[UoCt+Ti.t)]"’/ =U.(b,t4s)x 

Detailed demonstrations have been given in the general case / 1 6/ 


In the case of equation Tf, ^ one can define on one 
hand the propagator describing the exact motion of the 

particles in the fluctuating field £,* or U 

0(t,Vo) - L(£»t) Jj.-Ho 

( U operates on both X and U) 

or, on the other hand, the Weinstock progagator associated with 
the homogeneous part of equation (E.9A) 

\ A(t,U) = [ Lit) +'fcL’(t)] A aw si.M 

J 
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Here we have to stress the -fact that the fluctuation operator 3> 
equally acts on IJ and on A ■ The-solution of equations 
t) can thus be written 

ft — 

f 'ftj* A (t,U) + dJ:' Alt.tJ ]T.^«,a 


U(t,M f 7 bo) + 


U(b,t’j L'ft')F(fc‘j - U(t,b') F ^t')> 


T^Zh. 


wK«f« L* Lb*) =-E ^ i i 


Genera l ^igg£l --^Sag :£ii g5 

The kinetic limit of long times can be justified 
by the subdynamics method /13,14/; it consists in neglecting 
the influence of the initial fluctuation f /to) and in taking 
the asymptotic long time limit to — ^ . In this case we have 

f'/t) : — > I ds L' (b-6) f ^b-5) 

' c 

which describes the general mechanism of creation of fluctuations 
from the average function, by means of the fluctuating fields. 

On equation (l-izt) this limit gives 
C(b)54^L'(tlF'(b)) = dZ ^L'/bj U(b,b.-t)L'a--&))pM - dc <L'/fcjU/b,b.T^))e(t.o’ 

Jo 

which remains an integral equation. By decomposing 

u = u + u’ 


we have 


dlt)= <L'(t)U(t,t.»)L'(U)>f{t-») + 

Jc 

- jdi <:;L'(t) u'a,t-»)> c(t-n) 


TT .A& , 


Since we have 


4Ut)U(t,t')>' +4L'itiu'/t,t')> 
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the last term also writes *.s 

c It) . ( < L ' m L' (t')> f iw 

iflO t- 

' / J**' 

By using the complete (but formal) solution(E-*2A^ 
in •y,$ we obtains the Weinstock master equation 

+ ldt'^L'ii:}A(h,k*jV(e))f(Vj +4L’/tM(t,U)p'/to)> 


the kinetic limit of which gives the general kinetic equation 


^t^^t) = L(k)p(b) + 3^ (L'(k) Aa,k->) ^ I 


in an explicitly non— Markovian form d^^ends onf(b-?«) at 

earlier times) . The second term in the r , h , s . represents the 
turbulent collision operator in its most general form. 


The subdynamics method allows one to directly obtain 
the same equation in an apparently Markovian form { including 
however all the non— Markovian effects : see A6 /) : 

lU s + [dis V If (t) 

= [C(t| + &aij f tt) 

These two forms are rigorously equivalent since the turbulent 
propagator V is defined as the propagator of the kinetic 
equation itself: 

V It.U) S 1^ L (b) + ^fbj] V(t;bo) T.22. 

thus — — — 

f f ,U) 

Equation I!' 22 remains a non-linear equation for the turbulent 
propagator. This dynamical non-linearity can be expressed also 
in the non-linear equation for the turbulent collision operator 


i 


zed 


3-10 




At, <L'(t) V (t-i.t) 




where V is a superoperator/ non-linear in the operator G : 

V -- X+ (.'"ik' Lk(t') + Sit')] r.?r 

This kind of operators have been studied in detail in Ref ./15/. 
The formal solution of(%-tl) can be written in terms of the 
free propagator in the average field 


'It o'a.k.) = L(t) u'ii.to) 




in the form of the DYSON equation /1 8/ 

V(t,k<.)= Git) V(l'.t-) 


i. 



where (y is -jiven by JT.ZV in terms of A and V • The ranonrall- 
quasi —linear approximati^^QjJ^hich is equivalent to the Direct 
Interaction Approximation (DIA) introduced by Kraichnan, 

consists in approximating the Weinstock propagator /\ by the 
average turbulent propagator V ■ One obtains in this way: 


K9U 


.it' 



ir.zi 


where the RQL turbulent collision operator is defined by the 
following non-linear equation in terms of - 



]T.i3 


Equation ^E-2o,Z») describes the kinetic evolution of 
the distribution function ^ s -3. priori it involves all 

collisional effects. It has however the same apparent form as 
the equation obtained usually by introducing fluctuations 
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in the collisionless Vlasov equation , although such a procedure 
is not clearly justified usually. The equation obtained here is 
the actual justification of this usual procedure. Here the effect 
of individual particle collisions only appears later if one takes 
into account the specificity of the Klimentovich equation ; this 
appears in the binary correlations 

f.x= < 

which involves the "self correlations" 

r TC.io 

in addition to the usual ("distinct") correlations/ 19/ : 

p' p Set-F* 

For pratical purpose, it is sufficient to neglect self-correlations 
in order to describe the so-called "collisionless plasmas" from 
the Kllmontovich equation. Correlations are then turbulent 
correlations, and the equations are identical to those obtained 
by the common procedure of the "fluctuating Vlasov equation"but 
the justification is more clearly exhibited. 


II _^ E._ Quas i- linea r_ eg-uatio n fo r '* col lis l on les s ” tu rbulent p lasmas 


Approximating the general kinetic equation TT, 5 o or 
H consists in approximating the propagators. The quasi-linear 
approximation consists in retaining in the turbulent collisions J£ 
only free trajectories ( in the average field) determined by X.Z6 
One obtains 




3.32 


where 








flpo 

d5, < L'(tJ 


\0 


The free propagator U is given by 


j 
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^ Lit) U“a.U) T-3H 

i.e. ^ 

UVt.tc) , «t= j d*' l-v-2 - < JT. jr 

to 

and, in the absence of any average field ( I.e, for an homogeneous 
ensemble) it reduces to 

, -n-to)Y.y 

Uc^-t|to) s e. ir-36 


It is important to remark in the collision term It.'jj 
that the profs^ator takes into account ( here in an 

approximative way) the non-Markovian feature of the general 
kinetic equation. The effect of the propagator in the 

electric field correlation consists to introduce a Lagrangian 
correlation, here taken along a free motion i in a renormalized 
theory: alon the average motion) . 


The explicit form of the quasi— linear equation can 
be easily obtained in the absence of an average field. P/e then 
have 

L ^ - V .7 

(j*Ct|bc) —9 U<jlb|teW ^ 

and 




— — 

■^o 


+ 5 v.y - 


nr.^? 


Due to the obvious non-cotx'jnutation 




^ V.V 


e 



Tr.>s 


a gradient term appears (beside the usual diffusion term in 
velocity space ) which comes from the nori-Ma'rkoyian feature of 
the kinetic equation ; 
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fmuesBCStasKsa* 




where 


Pcv,t) = 


fp. 


ai |'Cx.'»y,t->))= 




^Ic.V? L 

dif C S (-5) 

- =tc 


TT.<lo 


and 


FtY.t).- 


d-6 S <E'^Y.tl E'Cx-^v,U)). 


€ll ■& 

o 


dk" 


-Clcv ^ 


U| 


- Vc 


jr.w 


This last non— Markovian term has been shown to b e of importance 
in presence of a strong magnetic field /1 9/ . Here the fluctuation 
spectrum is given by the Eulerian correlation ( assumed homogeneous 
here for simplicity) : 
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III. S0LnSIQN_QF_THE_ASYM£T0TIC_QUASI-LINE M-E2 nATI0N ; 
the_propagator_method 


The solution of equation II . is written in terms 
of the turbulent propagator V® 


f V*(k,t.) f Cs.y.t.l 


which is defined by the same quasi— linear equation: 


TT.^ 


+l.F(v,t).vJvYl.W iTi- 


where 5 
* 




. The formal solution is a 
time-ordered exponential operator: 

b ^ 


V '(t.U) = a; Ut' [-y.y tY.PMYY 


|TT.3. 


to 


/16/ 


where X^. is the operator of time— ordering which prescribes the 
ordering of the different factors coming from the expansion of 
the exponential , in the order of decreasing times from left to 
right. Non-commutation of the different terms in the exponent 
can be taken into account by means of an interaction representation. 
In a general way one can show /l 6 / that the solution 


U(t.U) , 1^' Ifttf) i-’fcavj 

•J I© 

of the equation 

(where A 3nd 3 cio not commute) is given in terms of the 
(assumed elementary) propagator 

9 Uftit.u) - Aai 

iby the following formula ( see aq. 3. 10 in Ref./iQ/) : 


H.fi 


•g;.y 


3 
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UJt.t.) X,«t- 


ft 




u 


By choosing As-'Ji.? we obtain 

^ L_ 


which can be put into the form 


_ -SY.V 


^ ^ X5. 


or 


-sv.v *■ 

V 5 e ■€- 


The general form of the coefficients 1 ^ ^ ^ has been 

given in eqs (5.14—18) in Ref. /I /. In the simple case of 
a stationnary spectrum, V and ^ do not depend on time. If 
we neglect in the exponent terms quadratic in the spectrum ^ 

( i , a , the weak coupling approximation) the result can be written 
in the following sisnp is and tractable form: 


a7”/ll ^ 

V It.Ujce e e e ^ e * 


where 




! "5 




■nr • U 


^ . 13 


-fCD.jWV"-') *■-4’'’ 


irr 


S., U,s) . ^ 


i« . /5“ 
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The physical meaning of these various coefficients 

obtained in the following w^. By integrating 
the exact equations of motion 

I Alt)- vft) 

jvffcU E[x(t|,t] .ic 

up to the order E^<v ^ and by taking the long time limit, we 
obtain : 

<X r X - 5 V 

■3 3 'BVx 


viCt-W> = V 

q: . >9 


H - 2.0 




We see that represents the turbulent contribution to the 

average displacement , the velocity, dependence of the coefficient 
^ represents an average acceleration or slowing down : these 
two terms are due to the velocity dependence of the diffusion 
coefficient. The 'term represents the elementary process 
of the quasi— linear description: diffusion in velocity space . 

The term % represents the associated spatial diffusion ( 
the Dupree damping term) and the term J represents 
an often neglected correlation between position and velocity. 

This whole set of terms have been discussed for the first time in 
Ref. /^ /. 


■m . 1 ? 
Ta • 



3-17 


Inversely on can check that the propagator obtained 
in '^,AA allows one to recover results , coming from 

the equations of motion, by letting the obtained propagator acting 
on X rX(lrj owl VsV^t) for instance: 

The obtained propagator thus takes correctly into account 
( and allows to describe) turbulent trajectories in the quasi- 
linear, approximation , t.e. linearly in the spectrum. 


In this weak-coupling approximation and for a 
stationnary spectrum, we obtain the solution of the quasi— linear 
equation in terms of the initial value as 


= e 


V /o) f C X ,V,o) 


■0 . 22 . 


Let us introduce the Fourier transform by 
— f . -Cir.X •— 

f Ut e ■ ■ 

f (V) . 1 (dX £ ^<'‘■''1 

'i En‘J> J 


We thus have 









m zr 


or 



nr. ^6 


3 
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transform 


For later use, let us also Introduce the Laplace 

f -C to t _ 

f Tto’) c \dJr e f (fcj 


We obtain 


_ fpc UlU).'k.Y)Z. -C|c.jV(y,5) -lc.|{y,?l.k \M(m\ 

d% -c e e e e 9 <v,Vs«) 


nr.2? 


In the simple case of a stationnary spectrum the 
solution of "5.2, can be written 

— , -tv.V + t .[^5(y)rb +F^y).Yj 

V*Ck,o^ = 6 ."’'^'"’^ " ^ 


and the solution '^•27 can also he written formally as 

f^Cy.io) r die e ^ 

' o 

KLoh -vU.y S + - 

= e e 

■'o 

that will he compared later to the non-asymptotic case. 


3 
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IV. RENORMALIZED QUASI-LINEAR KINETIC EQUATION 


Ranormalizing the quasl-llnear equation essentially 
consists In taking Into account the turbylant propagator In the 
collision term. Instead of the free propagator. One thus have 
( see 3.24 ) ! 


with 


Gtb). 


o. 




<LWfc) L' (!•-&)) 


■ ff .2 


and 


V(t.WU U"Ct,to)^- 


t 

dt' 

io 


iT.3. 


This last equation Is the PQL or DIA approximation of the DYSON 
equation ■ It Is equivalent to approximate the Welnstock 

propagator A Jby the turbulent propagator V . In equation 
the present approximation Is equivalent to keeping only 
the first term In the r.h.s,, (next term are thus responsible for 
"clumps" and other higher non— linear effects) . 


The weak— coupling approximation of eq. q.2 for the 
propagator consists In considering in_this equation only the 
quasl-llnear approximation of the collision term: 




< L'ft'; L’ct'.fi) v’c t-s, t'j v 1 1', w 

0 Z.9 


This Is the 30URRET approximation of the DYSON equation /II / . 
It Is equivalent to 

and the solution 


j 
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-'to 

has been developed explicitly in”wI.lo-H , 

The important physical effects introduced t>y this 
renormalized approximation is the main motivation of the present 
comparison between the techniques of propagator and Green 's functions . 


j 
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V. NON-ASYMPTOTIC QUASI-LINEAR EQUATION 


The asynptotic form of the general kinetic equation 
or of the quasi— linear one "it has been obtained either 
by the projector technique and the subdynamics method, or by the 
usual arguments: 

l^c) o 

m H 




Seme authors however consider the non— asynptotic equation obtained 
by sinply neglecting the initial fluctuation in the general 
master equation. In spite of the fact that there do not exist, 
to the best of our knowledge, any rigorous justification of the 
consistency of such an equation, one can consider that the long 
'time limit may be taken afterwards; such an equation presents 
some advantage namely for introducing Laplace transforms . 


In the master equation we simply take 

Vo c o 
» o 

and limit ourselves to the quasi— linear approximation 

Ad.t'l ^ u'Ct.t'J 


1 




X- 3. 


We then obtain the non-asynptotic quasi-linear equation ( 0) =0) 

. d'5 ^ ,t) e'Cx.TiV,h-'&)^ < ^ C i 






Jo 


where the integral kernel can be written in terms of the (Eulerian) 
spectrtim : 

, t --sy.? ^ 


\< (t): t . at € 


S e 

'I 


<T>V 


We assume here an homogeneous spectrum ( S is x-independent) . 


i 
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In Fourier tcans form { 1 h,^) ) we have 



- a'W.jc 


e 


t 




d5^ • e 


'i‘L.V'6 


I- ^ 

“i £.fe 


or 


Q)^+Ai£.v] ^ 


ft 




<^v 

- Jo 


5I . d-6 HCv.^-.s).^ 

K. /-v\/ s: ^mv 


where 


Kh)-.^- . 

q>V V 


■€. S 

=£ «T>V 


= . H ( W,v,-5) 

ov = q>v 




v.y. 


The so-obtained equation (V,? ) exhibits both 

non-Markovlan ( |l (t->) 4 "f(tj ) and non-asvmptotlc ( ) 

features. By neglecting these two effects, the turbulent collision 
term would reduce to 


-Ui h 




<i>V 


v-9. 


where 


'fin 

Cl'S H(kpo,v,^), 


o 


€t? €. " " = D IVib) 

- -£ = 


\)Mo. 


gives the usual velocity— diffusion term In the homogeneous case 

{ b = 0) . 


Let us Introduce the Laplace transform of the 
complete equation { V.9 ) : 

9 (u.)= 


t 


and let us restrict ourselves to the case of a time-independent 
spectrum. Equation v,? then writes: 


4 



V.IZ. 
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1 j^- A to 


+<u,i]f^(5i.w) = k-/u))f|,(v,u)4p^(!f.u.) ; 


Bo»3eBwe«Me»»«es!»B«*»»'*Bfs*"®s»''^ 


where 


k(u.) r 


k(*) .!• H 




\ 


and 




>» A:{w,is.v)5f yCivS 
d’s -e ^ 


S (it) = H j)i) 


such that 


H (u.-lsY ■511;) 




JJV 




3 
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VI. SOLUTION OF THE NON-ASYMPTOTIC QUASI-LINEAR EQUATION; 


GREEN'S FUNCTION 


VI.Aj. _ The Green_f un ctio n given b g H ORTON 

By writing equation n in the form 


- A J_U),W.v -A r 




tL.-i,. 


one can define the associated Green function by the equation 

-A r Lo.lcV -A k'/toH Q -- S^V-y') W.L. 

U>M 

and the solution of V./t6 can be written 

^ f ^ ? 

' ^ t.-o) 

- tc,jr - 

The solution of equation yj , z has been given by 
HORTON et al. / F / in an approximation which is equivalent to 
the weak— coupling renormalized approximation used in Section jlT , 
in the propagator formalism for the asymptotic equation. Our 
aim is to compare the approximations used in the two formalisms. 
Their solution writes 4 ui 'Ref •/5'/J •• 


G i\LX')- 

io,lc 


r^> 


— KLtL?^ A Tie - i ^ ic.i/ 4 w.v')is * ( v-w'j^ 
:! c . p ViF 


This form actually assumes an isotropic tensor 
V. H .V =■ H V*- 

W . H .W c H 


ii.r 


aH 


This is equivalent to neglecting the velocity dependence -L' =0 

«^v 

as well for the dependence of H in lo-k.y as for its dependence 
iny . 


3 
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In the framework of these approximations , we will compare in 
Section^ witir'the asymptotic solution obtained from the 
propagator 8.26. 


Let us introduce a Fourier transform in velocity space: 



zn^ J ' 


■^.6. 


The solution Vl.'i can thus be written in the considered approxi- 
mation ; 


fPO 


f (v,to)s dv' di — 

, r 


~ L 2 -- - J 

e e 




<4 « ' 


'ai- 


The dil' integration can be performed (v*rV-W) 


Av' e 






1 2 


dw €- 




, „ -SiCvik>>) 

(^itnU-c) e e 


'T7?' .y 


by using 


^.W _\>W^ 


idw ^ 


- 


gV'il’ 


W.JI 


The solution V)4 




can thus be written 

xujs --ik.yt ^y. 


dt €. 






e 


4 
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By collecting terms we finally have: 





X(\y ACOS --Clc-vs, 

_ - - - ^ 

€ e e e. e -e. 


KJ.AA 


ximxmmms9atssnais:iftf^99t 


One recognize here, in a scalar form, the spatial diffusion 
term W ^ H (Dupree damping) , the velocity diffusion term 

{ H ) and also the crossed term 0| H 1< corresponding 
to the correlations between position and velocity. 


VI . B T he_Erop aq ator f ormalism 


Solution of eq. Ji A can be written 

l-'t 


f = 1^- A W.v -A 


A L^- ^ - - 


^ a" X 

-aW.vx 4 5.1 .VlfuA.y , v). 1^1 

dx e ' 

e 

e 

1 J 




Here a propagator appears 

- , ..-Ij.si vxl-H/uxtv,v).l 

U . (^,<») - e 


OM •= 


OV 


loA.v 






which is solution of 


\ 



1 Hfo.t.U,V>l l 0 (,,c) 

Uj_fcy 


and it generalizes to the non-asymptotic case the previous 
propagator ( for stationnary spectrum) : 


vs- If. 


j 
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- A W.v z +6 — . T 3) ly) P <y' • Y 1 

ov L-^- ov * J 


vj ff 


It can be seen that in the non— asymptotic case ( vi.'M{ ) the 
non— Markovian contributions are taken into account in a much 
simpler way as compared to the asymptotic case ( F terms in ) 
This is due to the tensor y/u>-kv,v) : its approximate form for 
to ® If.y is equivalent to— the Markovian approximation ( F *0) 
of the asymptotic case because 


The knowledge of the expansion of the non- 
commuting^ operators in the exponent of ( Vl.U ) : 

V U,o). e. e e - p - - - g 

7=ck 


which for F '= 0 gives 


VI. rg 


\ I ^ 

Veik 




VI It 


we iramediatly obtain here in the same approximation: 


U (^*0) ' e 

Uj.Vc.u 


■e, e 


The solution Vl.l2. can thus be written as 

; ^ ^ » d J J I 










VI. i« 


VI. Zt 
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VI. C. Comparison 


When compared with the expression VhA'i obtained 
from the Green function, this last result has the advantage of 
taking into account the dependence of H upon both Uj-W.v 
and V ■ The Fourier transform in velocity space, which 
has been introduced to integrate the Green function actually 
needed to qeglect such a dependence. For‘ this reason, one 
can say that the simple result obtained from the Green function 
from the non-asymptotic equation is equivalent to the approxi- 
mation 

performed in the result of the asymptotic equation if we futher 
neglect the non— Markovian contributions (i.e, Eq. where 

^ ^ y 5 are given by . 12-15 with F =0) . 


In conclusion, we have shown that the approximation 
used by HORTON et al. to calculate the Green function Vl v in 
the case of the non— asymptotic quasi-linear equation is 
actually the same physical weak-coupling approximation used in 
the propagator formalism. However, the use of Green 's functions 
seems to be more delicate if we would like to take into account 
the tensorial feature of H 

“ jbut mainly its velocity-dependence { and thus the non- 
Markovian and non-asymptotic features) . 

Neglecting these two effects is equivalent to. only 
consider the asymptotic and Markovian quasi-linear equation. On the 
other hand/ the propagator technique allows one to take into account 

the tensorial feature of V and 5 in the asymptotic case 
{ or of M in the non-asymptotic case) 

their velocity dependence (i.e. average displacements effects: 
seeJi. ) and non-asymptotic contributions in the case of the 
tensor H 


3 
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SECTION 4 


A NBA/ KINETIC DESCRIPTION FOR TURBULENT COLLISION INCLUDING MODE-COUPLING 
J. H MISGUICH and C. M. TCHEN 


ABSTRACT 


The micro-dynamical state of fluid turbulence is described by a 
hydrodynamical system. This is transformed into a master equation in a 
form analogous to the Vlasov equation for plasma turbulence. When the 
total distribution function is decomposed into a mean value and a 
fluctuation, the evolution of the mean distribution satisfies a transport 
equation ^called the kinetic equation, and contains a turbulent collision 
that represents the statistical effect of the turbulent fluctuations, 
while the evolution of the fluctuation will form a transport equation fo^ 
the turbulent collision. Themechanism of this collision is investigated. 

The report compares our theory with the clump theory of Dupree, the 
renormaliration theory of Misguich and Balescu and the direct interaction 
theory of Kraichnan. 

The paper is completed by Dr. Tchen at the City College of Maw York 
in collaboration with Dr. Misguish at the Centre d'Etudes Nucleaires, 
Fontenay-aux- Roses, France. 

The work is the result of cooperation befween Dr. Tchen at the City 
College of New York and Dr. Misguich at the Centre d'Etudes Nucleaires, 
Fontenay-aux- Roses, France. 
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INTRODUCTION 


In plasma and fluid turbulence, the main average 
effect of the fluctuations consists to enhance (cr replace) the 
role of particle collisions and to enhance consequently the 
transport processes ( turbulent diffusion, turbulent viscosity. 
The origin of such phenomena is generally referred 
to as "turbulent collisions" . 

The starting point of theoretical descriptions are 
the Navier— Stokes equation for the macroscopic description of 
fluids, and the Klimentovich equation for the microscopic 
description of plasmas. The latter reduces to the Vlasov equa- 
tion in the case of "collisionless" plasmas, which behave like 
an incompressible fluid in phase space. 

The most widely used among the theoretical descrip- 
tions of such processes are the Direct Interaction Approximation 
(Dlk) introduced by Kraichnan in fluid turbulence, and the 
Renormalized Quasi— Linear (RQL) approximation introduced!^ 
Dupree and Weinstock in plasma turbulence. Both approximations 
have been shown to be ana,loguous . /!/ 

Similar formalisms can thus be used to describe 
both kind of turbulence, and we will adopt here the language 
of plasma physics. The results can immediately be translated 
in the fluid language by replacing electric field fluctuations 
by pressure gradient fluctuations 

Turbulent collisions effects are described in RQL 
or DIA approximations by means of a Lagrangian correlation of 
fluctuating fields between two points whic'g are separated in 
space and in time. These distances in space and in time are 
are actually related to each other by a complicated trajectory 
which involves all the dynamical problem; in RQL or DIA this 
trajectory reduces to an average diffusive trajectory. 
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The purpose of the present work is to analyse the 
effects of higher mode coupling terms on this approximate 
description. Dupree for instance has introduced so-called 
'beta terms" in plasma turbulence, the role of which is to 
ensure energy conservation to a higher degree of precision. 

Such terms have been proved to be important in drift wave 
turbulence /2/. 

The introduction of such higher order corrections 
has however been made in a rather intuitive manner ; a simpli- 
fied approach is presented in Section 11. Here we derive new 
results for turbulent collisions by using a general kinetic 
formulation of plasma turbulence (Section III) , Moreover we 
prove that these higher order corrections already appears 
from the non-linear dynamics, even when the non-linearity 
introduced by Poisson equation is not yet taken into account : 
they cannot thus be reduced to self-consistency effects only. 

Our main result (Eq. 111.14) consists to describe 
the deviation of the turbulent collision term from its RQL 
or DIA approximation in an exact way. This deviation involves 
a generalisation of the Beta terms introduced by Dupree, which 
includes here all higher order mode-coupling terms. A new 
approximation is proposed in Eq. (111.20) which goes beyond 
the RQL description; an approximate treatment of this equation 
is left for a future work. 

Another approach is presented in Section IV, which 
uses a simple non kinetic treatment of the Vlasov equation. 
Although the resulting expression involves four main contributions 
(Eq. IV. 24) , we have been able to prove that in the kinetic 
regime an unexpected cancellation occurs between the last two 
contributions, and we recover the general result of the kinetic 
formulation . 


The consequences of this cancellation are examined 
in Section V where we focus our attention on a quantity 
defined as the infinite time integral of the Lagrangian 
autocorrelation of fluctuating fields. Exactly like in the 
classical case of Brownian motion, we demonstrate here that 


j 
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tn an asymptotic description, goes to zero for an homo- 
geneous and stationnary turbulence. Of course, when different 
scales can be introduced to describe the fluctuations, the 
Lagrangian correlation of the small-scale fluctuations remain 
inhomogeneous and non— stationnary, and this avoids the 
corresponding to vanish in the asymptotic case. In such 
a scaling description of turbulence, C£ can be used as a bare 
description of the turbulent collision : higher order mode- 
coupling terms then will introduce a dynamical shielding effect. 
This description will be developed in a subsequent work where 
an explicit calculation of the intrinsic or shielded turbulent 
collision, will be performed. 


j 
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II. TURBULENT COLLISIONS AND HIGHER ORDER TERMS 


II.. JEkaheaiaJtJ.aal. Jjit£oduc±ion 

The recipe for the introduction of the higher order 
Beta terms can be summarized as follows. The Klimontovich or 
Vlasov equation for the distribution function f(x,v,t> 

\ ^ f 

with 

L c-v.v - 1 erx.vvl 

„ 3 

( V = ^ , q is the charge and m the mass of the particles, 

E. in plasma is the electrostatic fluctuating field) can be 
separated as usual in fluctuating and average parts: 

■ -it f'-. If' + Df + SL'f 

with the usual notation 

f ? 4F>-- R F 

and B=l-A denotes the fluctuating part of everything to its 
right . 


The free propagator in the average field L is 

defined by 

UJt.V.) = Lit) U.tl.V.) 

11.4 


In the simple case of linear trajectories this reduces to 

and the solution is a simple finite displacement operator in 
x-space: 

= e 

11.6 


j 
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or 


Solution .of (II. 3) can be obtained as 

u^ct.vj [I'a'jfa'j 


II,7a 


f ' . U, [L'F +3 l'F'] 

II. 7b 

where time dependence is taken into account by the notation i(r 
which denotes the time convolution. This solution is formal 
since it is an integral equation which remains to be solved. 

Substitution of (11.7) into (II. 3) is equivalent 
to an iteration and gives 

' II. 3b 

A simple handwaving argument can be applied in the last term 

in the r.h.s. : in '5L'P’=L'F'-4L'F'> ' the projector 

B means that the average has to be substracted from L' f ^ • In 
the final calculation of average quantities the L* cannot 
be taken in average with : simple averages involving this 

last U can only be taken: 

i) either with the other U which yields 

ii) either with other fluctuating quantities which will appear 
to the right. This case can be schematized by an arrow 

4l' U.»^L' F’> 

^ — ... J 

which indicatesthat L* is operated by Uo -but is excluded 
from the bracket average. 

The last term in (II.3b> thus involves at least the above two 
terms ^ and (II.3b> can be written: 

\F'=LF' + i-'f t<L'y,*L'>F' + <CL'Uo*L’f'> 


Both underlined terms can be used to define an average renorma- 
lized propagator and the solution 'writes. 

I^L'F + •BL'U^X-L'F -H < L' L 


II. 9 



J 
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When substituted into (II. 2) this result allows us to write 
the average equation __ 

V 4 L' l' Uo L' f ’> 

t — ^ 11.10 

In the r.h.s. the first term describes free motion in the 
average field, the second term a turbulent collision term 
(here in the weak-cdiupling RQL approximation) . The third and 
4th terms are higher order corrections to this approximation 
of the turbulent collisions. The complete expression for the 
collision appears to be in this handwaving derivation: 

<L'VVU'>f 4-4L'V‘’jfBL'L),ifU'>F 

+ <L'V%L' U. 

^ ^ 

11.11 

Actually the last term ( equivalent to the Dupree Beta term /2/) 
is of the type 

<L'L') <L'f> 

which involves the complete <:(k) ; Eq. 11.11 still remains an 
integral equation for , which describe some kind of 

renormalization or shielding of the RQL— weak coupling appro- 
ximation given by 

c'’ . < 

The purpose of the present work is to elucidete the role of 
such higher— order implicit terms in the equation which defines 
the turbulent collisions . 

This above schematical introduction of the higher 
order terms can be made slightly more precise by considering 
the wave— vector dependence in Fourier transform. 
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For the case of a turbulence which is homogeneous 
In average, we have a space-independent average function 


and 

and 


f iY,t) 
L 5- V.5 


AW. —K «DV 

since we have no average Vlasov field. In this case (II. 3a) 
writes 


11.12 


In order to perform the iteration, we write the formal 
solution for ■ Neglecting as usual the initial value 

term, the asymptotic solution writes 

t 


' K.k' 




cU:' e 


] 




11.13 


and (11.12) can be written in the form: 


4k' bU' (» 






11.14 


Let us consider the structure of this general expression. In 
the last term, the first B to the left prescribes the fluctu- 
ating part o*f what follows. (The integral can be replaced 

by the discrete summation with appropriate factor 

which are ommited here for clarity) . In a homogeneous system 
this implies k" / 0 .In the same way the second B implies 
W+W ; two values of 1^*' appear to play a special role: 
i) k-"*W which gives L*j,(k*) 

which is analoguous to the Dupree Beta term 
ii>k“**-W‘ which gives 

i.e. a phase coherent term. 


j 
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Neglecting higher non-coherent mode-coupling terms gives: 



11.15 


Here the k' summation is an approximate form for the complete 
mode- coupling term "bL'f in II. 3 a or 11.12 •. 



11.16 


Let us consider the structure of this equation: 

” The term /‘■V appears as a modification of the source 

term in P ; it has been interpreted by Dupree as a modifica- 
tion of the average function p into an "effective" average 
function ! 

r.f 

” The term appears as the weak-coupling approximation 

of the turbulent collision term ( coherent part). 


In Ref. /3/ these both terms appear as 


Tbur-^C.f, . Ct - ]dt' L4f) C;u-t') 


-t 


4> 


11.17 

i.e. a phase coherent term ( turbulent collisions) and 

the Beta term ( ). This corresponds here to the turbulent 

collisions 

- - te* - - - 

and the Beta term 

rv 


11.18 




li’ L' IV) C^(h.h') -- dt' 1 _ L < i' lil 

5 ^ ^ )r • ~ 

Jii - S 


The role of the term L' . L' r 

I Vr-k* r 

not apparent . - - i 


11.19 

in Dupree's description is 


a 
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III, EXACT KINETIC RESULTS 


The appearance of turbulent collisions and higher 
order terms can actually be described more rigourously by 
using general expressions which have been obtained in the 
framework of our kinetic theory for plasma turbulence. 

III. A . General kinetic equation for turbulent p taamaa 


In the kinetic regime of long times the exact 
solution of II. 3a for the fluctuations can be written /4/ 



This is a basic formula which describes the fluctuations in 
terms of the average distribution function. It holds in the 
kinetic regime of times long compared to the correlation time 
of electric field fluctuations ; in such regime the influence 
of the ••nitial- fluctuation has been shown to vanish exactly /4/. 
Here A is the Weinstock propagator defined by the homogeneous 
part of Eq.( II. 3a) : 


: [udl +-6L'^t)J A^t.Vc) 
with AU.,VoV- AU.t) 


111. 2 

111.3 

111.4 


and where is the time— ordering operator /5/ which prescribes 
that in the expansion of the exponential operator in power of 
its exponent, the various operators have to be time— ordered in 
order of decreasing times when t t 


J 
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Solution 111.1 can be checked immediately: by using (III. 2), 
Eq. (111.1) gives (II. 3a) (see 111.4) . 


This simple exact solution allows us to write 
the average equation (II. 2) as 


f = L Utl + 6 (t)] f (t) 


III. 5 


iSCUt 






equation At self: | 

111.7 

Vdst.)-- e. 

111.8 


This operator allows us to write the exact solution as 

f V (t,f) 

III. 9 


i S 5 rB' r -i fr ei»a fc ± on 


Let us now perform a kind of iteration of Eq. (III.l) . 
A simple integral relation can be obtained between the Weinstock 
propagator /\ (III. 2) and the turbulent propagator (III. 7) : 


A V + V a - V t a =» ^ 


III. 10a 


j 
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i.e. 


(tw 




III. 10b 


This integral equation can be deduced from (111.2) by 
treating t^L'- G) as a perturbation: 


By substituting (111.10) into (III.l) we immedia- 
tely obtain the following iterated formula for the fluctuating 
distribution function: 

(Oo «. 


dl% 


V f tt-6) 

B 


B "B 


111.11 

By substituting this exact solution in (II. 2) we find a new 
expression for the turbulent collision term: 


CU)s<^L'F'>= 6ttl^it)-. 


dT> U'tt) F^ty 


111.12 


which is given by 


at ) = 


9» 

dz 

o J 


% 

dJ9 


< L! tt) V (t ,V-i) L' i V-v)) F l-Ki) 

/( L'(l) V UiV-») L^L'(te) -&{t-60 L'W->)) 


111.13 


3 
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Introducing the variable Ss S-& we have 

do Cdo (9o 

fit's S \<i.e 
Q Jo "'o 


and finally the result can be written as (0 S ) : 



111.14 


The first term is nothing else than the RQL approximation of 
the collision term, which has been developed in plasma turbulence 
in analogy with the DIA approximation in fluid turbulence/6/ .This 
approximation is obtained by replacing in the exact expression 
(111.6) the fluctuating trajectories described by the Weinstock 
propagator A by the average turbulent trajectories described 
by the turbulent propagator ^ . The second term in (111,14) 

takes into account the deviation from the RQL description, by 
means of 6-) 

1 1 1 J.C I. ^H ig her _ordgr_ terms ^Jbeyond-RQfe 

In order to make a bridge with the more intuitive 
calculation of Section 11, it is interesting to put this 
correction in another form. From (III.l) we have indeed: 

■'o 

111.15 

Since 

f (.V-t-O? V f rV(V--5-S,t) 

111.16 

we can write (111,14) as : 



J 
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C(tl= |ck) = d-i <L'(t) v/t,M L7v-ii)V^u,t^ p(t; 




- fdx ■^L'lt) Vit,V-») S-tt-t) 


111.17 


When ccxipared with 11.11 the present exact result (111.17) can 
be analyzed as follows. The first term is the RQL approximation 
of the collision term, while the first term in 11.11 only 
represents the quasi-linear result ( A U« ) . The second 

term in (111 .17) is a generalization of the second term of ( 1 1 . 1 1 ) 
i , s . a term involving at least three fluctuating fields; here 
A remains a fluctuating quantity which may introduce corre- 
lations between four fields, a.s.o, fhs third term in (III. 17) 
is the general form of the Beta term obtained in ( 1 1 . 1 1 ) . The 
present general formulation avoids using approximate treatments 
of averages like that denoted by the arrow in Sq. (11,11) and 
avoids elimination of higher order mode— coupling effects. 

Moreover it allows to exhibit the implicit feature of this 
higher order Beta term: the deviation of G from the RQL appro- 
ximation of G- involves a Beta term which depends on G 
itself; this was already apparent in Eq. (11,11) where 4 L* p 
was given in terms of p')’ 

In summary, we have shown that the turbulent 
collision term involves a first lowest order approximation 
(QLin 11.11, RQL in 111,17) plus higher order terms which 
have been calculated. This latter correction involves an 
important , so-called Beta term, which actually depends on 
the complete collision term ( "implicit'^ terms ) . 

^Aggroxlmations 

Two levels of approximation can be used to 
describe (111 . 14) . 


i 
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1 ". 


The RQL approximation for is solution of the 

non-linear equation (replace /\ iby V in III. 12): 





Jo 


III. 18 


with 






III. 19 


2® . Higher order corrections can be obtained by replacing 
^ by V in the iterated formula obtained in (III. 14) . 



This equation is the main tractable result of the present 
work, It represents a highly non-linear equation for the 
collision term of the average equation 

f = Lf + Clt) 

described in an approximation which goes beyond the RQL 
description. The calculation of an approximate solution for 
this new closed equation i s left for future work. 


i 
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IV-. NON-ASYMPTOTIC TREATMENT FOR TURBULENT COLLISIONS 


An integirall equation similar to (111.17) 
for the collision term of Eg. (II. 2) has been obtained from 
the Vlasov equation directly / S/ . We will show here that the 
physical content of the different terms appearsto be quite 
different . 


iy« .Integral eguatlon 

The fluctuating equation (II. 3) can be written in 

the form 

* l'f 

IV. I 

This form naturally introduces the exact propagator U defined 
by 

'5tU(t,U= uO.V.i 

IV. 2 

uawx e 

IV. 3 

Contrary to the Weinstock propagator, the exact propagator 
(J does not prescribe all intermadiate states to be fluctu- 
ating ones: |J allows "transitions" from an average "state" to 
a fluctuating one, and vice-versa. The Weinstock propagator A 
(III. 3) in turns describesthe irreducible propagation of 
fluctuating states /4/. 


Solution of (IV. 1) can be written ( neglecting 
the initial fluctuation) : 


F' rU*L'f - Ui(-4L'F'> 

i.e. 

ua,vj [L'aH(f) . <LVf) FVfj>J 

A. 


IV. 5 
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IV. B. Three contributions to the collision term. 


Let us transform Eq.(IV.S) in order to introduce 
explicitly a first term which in the asymptotic limit (t^ — ) 
gives the RQL approximation. For this, we simply decompose 


U = U -i-u' 


IV. 7 


where the average U of the exact propagator actually tends 
to the average turbulent propagator (III. 8) in the limit: 

of time intervals long compared to the correlation time of 
electric field fluctuations. Then (IV. 6) becomes: 


ciu.c,jtn 


where the three contributions are given by 

I 



In order to build a bridge between the present 
treatment and the asymptotic results of Section III, we have 
to calculate in terms of A and V . For this, we have 


j 
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to derive from the definitions of the propagators some simple 
exact relations in the form of integral equations. 


Let us first derive the (integral) relation 
between U and ; 

U . U, - U. *■ L' U 

IV. 12 

which is readily obtained from (IV. 2) and (11.4). We have to 
note that in such relations, the two propagators in the second 
term can be equally inverted since both formulas give the same 
iterated explicit solution ; 

u = i (u,*lTu. 

IV. 13 


Let us now express in terms of V and 

From (III. 7) and (II. 4) we have: 

V r Uc + V Jf & U. 

IV. 14 

Combining this with (IV. 12) gives 

U=Uo+V »L'U -V )(-<SO.»L'U 

JV. 15 

Let us also express U in terms of A ■ Eq.(III.2) can be 
written: 

A t L + l' - AL*J A 

which gives 

A r u - U Pi l' A 

, IV.16 

Then IV. 15 gives 

U = L). +V^^L'A+ Vjf- -V jKiU.ifL'U 

JV. 17 

Here U© in the last term can be eliminated by 'using (IV. 12 
and 16) : 


Jl 
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U^JtL'u = U -Uo r A +U^fAL'A-Uo 

IV, 18 

and we finally obtain 

VJzU. +V j(-L'A t V JfU'U -V 

IV. 19a 

IJ. L' U if FM-'A -V/eG[A-Uo+U ;f Au'aJ 

IV. 19b 

The fluctuating part of U can thus be written as 

U'= VaL'A ifL'U ifAL’A-Vif'SA -VifGU jfAvJAj. 

l .e, 

U’ c V if -t- V U if A L'A 

Iv. 20a 


Or r explicitly 


le V L\:M A L 


d© va,U)[ibL’Ct-o)-(ia-eJ] 


ds U(t- 8 ,t- 0 A LYt.&J AO-S/t-&) 


B 


IV. 20b 


IV_j_D^_As]^gtotic_cancellation_of Ciw 


From the exact relation (IV. 20) we can write the 
second contribution to the collision term (iv.lO) in the form: 

j ^ L (t)(jj t U*3 L C t f t CjT2, 




IV. 21 


with 


u. 


Ctt tt), d^<L'(U 




IV. 22 
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KVo 

d% ^ L' (t) 


dsVft.^-e) [lbL'(V.p). 60 .ej] 


rz 

oisUfKeit'i) 




TV. 23 


The total collision term (IV. 8) thus writes: 




IV. 24 

Actually one can prove that, in the asymptotic limit, two 
last terms in (IV. 24) exactly cancel (see Appendix) ; 




JV. 25 


This implies that we have asymptotically 

ctt) — 


I Cjll) + Ct tti] 

bo P* 


IV. 26 


Now this result can be proved to be identical to 
Eq.(III.l4) obtained in the asymptotic description. This 
identity can be proved as follows. From (IV. 9) we have indeed 


U C 

too*«* 


uu 


fat <L'(t) va,v-t) L'a-M)fO-t) 


IV. 27 


which is indeed the first term of (III. 14) . On the other hand, 
we have frcm (IV. 22) : 



(PO {>» 


dB 


\di> 




JV. 28 


which is identical to the second term of (III. 14) . This achieves 
the demonstration of the equivalence of (IV. 26) and (111.14). 


We have thus seen in the general result (iv.24) 
that in the asymptotic limit of a distant initial condition, 
last twi terms exactly cancel. This means that in (IV. 6) the 
last implicit term is cancelled asymptotically by a part of 


3 



the contribution of in the first exPli^it term/ and we 
remain with the asymptotic result given by the kinetic 
formulation (III. 14). This remains a highly non-linear equa- 
tion for the collision term, even in the approximation (III. 20) 
where A has been replaced by ^ 
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V, ANALYSIS OF THE LAGRANGIAN AUTOCORRELATION OF FIELDS. 

V_j_Ai_La2ran2ian_autocorrelation_for_turbulent__glasmas__and 

§E2^SiSS-.S2£l:2S • 

Let us consider the exact relation (IV. 6 ) which 
has been derived for the turbulent collision term from the 
Vlasov or Xlimontovich equation: 



where the first term defined in terms of the Lagrangian 

correlation between two electric field fluctuations along the 
exact motion of the particle, which is described by the 
propagator U : 

V-Vo — 

c^tk) « = a-b ua.M 

AH^ *l»v/ ' “ “ ' 

- Jo 

V.2 

This autocorrelation of the electric field fluctuation can 
be visualized as the correlation between two points of the 
space— time (x,t) which are defined by the exact position of 
a particle which is at point x at time t and at a point 
x(V-'& ) =U(kiV*i))( at time t-t ( see Fig .l) . 



The time integral of this autocorrelation is fully similar 
to the well-known autocorrelation of fluctuating forces , which 





4-23 


appears in the classical description of Brownian motion. The 
friction coefficient of the Langevin equation 


r-1. F/t) 

dbfc /wi 1 - 

V. 3 

is related to the diffusion coefficient by the Eintein relation 


T 

and is given by the time integral of the autocorrelation of 
the total fluctuating force ( see p. 262 in /!/) : 




fit 

5VT Jo ~ 


V. 5 


Here the upper integration limit? is larger than the 
correlation ^of the fluctuations , but smaller than the 
relaxation time - 


S,<T 

* V.6 

The point we want to stress here is the following: the exact 
autocorrelation has a strong peak near the origine CS<5c) 
but also a long negative tall : 


(see Fig. 2 /8/) 


V.7 



3 
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which makes that the total intecrral ( from zero to infinity) 
of the autocorrelation is actually zero : 


P 


9 

The finite integral however 

r 

Jo 


v.a 


V. 9 


is a function, of r which reaciyja quasi— plateau value at 
(seep. 261 in /!/) 



VtB: — General- -resu-lt 


The point we want to show here is actually very 
similar. Let us consider Eq.(V.l) and calculate the second 
term. From the exact equation for U we have: 


U-. LU 

thus 

u = I u f 4 l'u'> 

This can also be written 

<L' tk) u’(t,ej> -- - L(t)J u itx) 

When substituted into (V.l) : 

CLt)- CJt) * - dk' CU'j = 


V. 10 


V. 11 


V. 12 


V. 13 


3 
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this gives an equation where fluctuating quantities have been 
eliminated: 

fh 






to 


V.14 


By using partial integration 

r t 


-lo 

we obtain 


tlV 1 Hft.l'jd it' Hll.k'J - H/t,k) 

Ot <T)h , 


V.15 


c^tt). ca) - fit' utt.t'j c(k'j 

“ kj 


V. 16 


By this way we have proved that 


res. ^ “1 

ft - , , . 


at’ Civ) 

to 



V.17 


This is an integral equation for ^ in terms of (Zi_ which 
is given by the time integral of the Lagrangian'autocorrelation 
of the fields. Here the integral operator has to be inverted 
in order to determine C. In terms of Ci. . An important- 
conclusion can nevertheless be reached for a particular 
homogeneous case. 

V.C. Asyrngtotic descrigtion_for_homo 2 eneous_and_stationnary 

In the asymptotic case ( tjj ) this equation 

4 L’(i) ua.M = h- 

m “ 


becomes : 

/■po 




urt,b-^)ca-Tv) 

As 


V. 18 


and 


As 


( dt' LL(t') + 
Jf.-6 


GCt’)] 


V. 19 


J 
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In the stationnary case t -independent, and 

CLU) becomes time— independent. We remain with 


c|;'aj = - Lit) 


dz \ICz) C 


For homogeneous systems, C and V are x— independent and 
L ’*-%V { average Vlasov field) . We thus reach the 

result that for an homogeneous and stationnary turbulence 
the infinite time integral of the Laaranaian autocorrelation 


of fluctuatina fields vanishes in the asymptotic description 





although the total turbulent collision term does not vanish. 
This result is very similar to Eq. (V.8) for Brownian motion. 

Yi2i_S2SSS3HSS2S-£2E-£HE^iiiS2l_22iii§i2S§ 


This result implies that 

Jt is interesting to combine this result with (IV. 25) 


V.22 




When reported in Eq.(IV.24) this gives 

+ Cj /t) 

^ — i i 

* CW - CU) +C(t) 


This means that in the asymptotic Analysis of the four 
appearing in (IV. 8) for the collision term ; 


V.24 

terms 





3“) .The total collision term is given by the result of the 
kinetic method developed in Section III ; 



V. 21 



These equation (V. 25^26, 28) are three independent relations. 
The importance of the last result is that it allows us to 
describe the turbulent collision term by a modification of 
the RQL approximation . The contribution of in (V . 28) 
represents the equivalent of the Dupree Beta terms, here in 
a general form. Their importance has been proved for energy 
conservation. We will come back later to this important point. 
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Let US prove the cancellation of Cnr ^nd in 
the asymptotic case. We calculate (IV. 11) 

s - I ckz ^ LVfe) (J 

'' A.l 

and use the exact asymptotic expression for C. (III. 12) : 

rT” 

■“ T--><»9 L 


Let us modify the variables: 


pT fT tr fz' . f&' 

at di' r cu' d-s r at' as 

Jo J-6 Jo Jo -^o Jq 3 

and eliminate ^ by the new variable Sr 6-s. We obtain: 

t r^' - 

Cj-(L - in' d4<L'/ui;'(i-,i-5'-*0)<L’0-5'«)A^f-tVs,t.i'jL'o-«';)Fif-5'j 

J 4 ) 

A. 4 

Let us now use (V.20b) for (J : 

(TJ'S f 

U'a,V-t'ts)= dV V(KV5)[liL'»-8)-6it-e)] i A^l••6,^-s’■*0 

Jo 

t'-S n 

V d.). Utu,k.>) AU'a-|.) ACk^il-sVs,) I 

•» - A.5 

The bracket can actually be related to the exact propagator 
U by means of (IV. 16): 

Jf.'&'d.c 


r A(t-e,k-t'4s)'i' d\> AL'{t-f) A{t-V,Kt'+s) 


4 



4-29 


("whereas b.t* ) . Substituting (A. 6) into (A, 5) gives 


u'aAsvo- 


[k'e Va,V-9)Il!.uVl-»)-G^‘-eQ 


A. 7 

When sqbstituted into (A. 4) this result allows us to express 

re Z 

^.S 


^ change into "S ) 


e|‘rn--<u;. - 

T-l* 


r 1 

rz ' 



VQ J 

o *' 


d'e < L'ft)V^t,V-sjLTbu’£KB).6^t-B)JuO-e.V-«<i)^ 

9 


A . 8 . 


On the other hand, (1^-23) can also be transformed 

by using 


Ss5.a' 

I i 

'0- Jo 

and we obtain ( :^S ) 


fZ 

eke 

o 


dl-6' r 


fitV 

o J 


ri5-'S' 

o 


A. 5 


Ci^lH . 

dz ^ LVt] 

•c ■ 

fits 


o •J 

0 


^ Va,Ke)['iLVV-e)-Si!ks)]ua-B,V-»-'*-r^ 


L' rKi-ti) AC^•*+4l'^-*'5L' 

A. 10 

By comparing with (A. 8) we have thus demontstrated the result 
(IV. 25) ; 

[as < u' (t)V 0:,v-9) 



r*-» 


C 

di; 

fits 

' 

o 

0 * 


4 L'(V-B-.S) A0-MS,v->^)L'a->.)> 


A. 11 
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ABSTRACT 


The hydrodynamical equations of turbulence are transformed into a master 
equation for the velocity distribution function. A group- scaling is introduced 
for the closure. The spectral balance for the velocity fluctuations of individual 
components shows that the scaled pressure- strain correlation and the cascade- 
transfer are two transport functions that play the most important roles. W 
derive this correlation and find a power spectrum for the horizontal components, 

while the. spectrum for the vertical component drops rapidly by going to the large 
scales. 
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In a strongly stable boundary layer, the spectral distributions of the 

_3 

large scale horizontal components take the power law k , while the vertical 

component decays rapidly toward the small k values. This spectral hump in the 

horizontal components are investigated here. We use the group-kinetic method 

1 

of closure developed by Tchen . Instead of the customary Fourier decomposition 
of a random function, we perform a group— scaling as a coarse— graining procedure . 

We replace the hydrodynamical equations of turbulence by a master equation for 
the velocity distribution, and exploit its advantage of being homogeneous and having 
Less nonlinearity, since the turbulent velocity is now an independent 
variable. 

11. BASIC EQUATIONS OF ATMOSPHERIC TURBULENCE 
A. Microdynamical State 

The microdynamical state of atmospheric turbulence is governed by the 
equation of continuity and the hydrodynamical equations of Navier-Stokes, as 
follows : 

7 k = 0 





(lb) 


A 

The total function, as denoted by { ) , consists of an average quantity, as 

denoted by 


= < } ■ 

and a fluctuation, or a deviation from the average, as denoted by { ) . !The 

/A 

evolution of the velocity lx is controlled by the kinematic viscosit 'y 3 and 


a 



a field 
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^ =1 £ 


( 2 ) 


which is the sum of the adverse pressure gradient 

A 

J P ^ 

divided by a constant reference density P , and the buoyancy 


(3a) 


£ =- 


fo 




(3b) 


The buoyancy is produced hy the variable density under the gravitational 
acceleration j ■ If the temperature T and the humidity Q have tliuir 
mean values T, Q, constant reference values , and fluctuations 

0 , , we can write, under the Bousslnesq hypothesis, 

f-fo ^ A f-T. (j>-Q. 

Q 




X 


(/.) 


The temperature and the humidity satisfy the following equations of evolution: 


-h - fr f = 0 

(\+ Q-V- v"']q = 0 , 

^ k>V 


(5a) 


(5b) 


with the molecular dif fusivityh) . 

The system of equations (la), (lb), (5a) and (5b) describes the microdynami 
cal state of stratified turbulence. 

B. Group-Scaling of Fluctuations 

For the study of the spectral structure, wc will bo interested in the Tlucina- 

'V' A 

tions U ,v , ^ 5 whose governing equations can be obtained by applying the 

(/ 

^ ^ ^ — . 

fluctuation operator Ab!~A to (1) and (5), where ^ is the 

averaging operator. These equations have terms proportional to the mean gradients 


J 



f 
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VT , VQ ^ 

treated as constant parameters. For the sake of simplicity, we combine the 

/\ 

tenperature and the humidity into a single variable 'V, called drift, such that 




The constant coefficient^ has the dimension of a frequency, and will always 
accompany the drift to form the buoyancy force. 

The differential equations which govern the fluctuating functions tC and 
can be obtained by applying the operator A to p.) and . These equations 
do not contain the mode-coupllng3, unless a Fourier transformation is made. 

However, the Fourier form, like the original equations for the fluctuations, 
contain too many minute details. A coarse-graining is necessary for a statistical 
study. For this reason, we introduce a group-scaling, by means of the scaling 

O 

operators A and A , which compose A, in the form: 

A = A% a', (8a) 

and the sequence continues by writing 

a' = ^ + A • (8b) 


Thus we have the macro-velocity fluctuation 

o A o 

A u =k u , 

the micro-velocity fluctuation 

/ /V / 

A u = u , 


the submacro-velocity fluctuation 

(Ij-v. (') 
A u = u , 


and the submicro— velocity fluctuation 


u = 


(9a) 


(9b) 


(9c) 


(9d) 
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O t 

The goups U. , U f (A, 'have their correlation times 

> t '" (10) 

ranged .in the decreasing order of magnitude, indicating their decreasing 
coherence. 

o 

By the aid of A , we transform (la), (lb) and (5) into: 

ulv.Tl u!,-ufr 

t & d ^ I' J Z'' 

7- iJi ^ 0^ 7- a^= o . 

^ *V 

It is not difficult to derive the following transport: equations of energy 
in the group form: 


(lla) 

(llb) 

(llc) 



with the transport functions: 



( 12 ) 


(13) 


(lAa) 


(14b) 


i 
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i 1 




\u 


=- i^yyy) \ 


(l/l<-) 


(14d) 


(14e) 


O J^O 

Thus we have the production functions i and iT, by wind shear V/ U , and 

o<W n/ ^ ^ ^ ^ 

Jty mean drift gradient ^ Ia^ , the buoyancy transport 5 the re-dis 

V ^ 


(0® " *T~® '7 

among the components <>( = 1, 2, 2>, th$ cascade transfers and I 

e ° c ° 

and the dissipations and c-^ 


tr ilm Lion 

— O 

Vc< 

The summation rule which applies to. 


the repeating Roman indices does not apply to the Greek indices. The transport 
equations have been approximated, by treating the mean gradients 

VlK , VW ( 15 ) 

C «V 

as constant parameters and by omitting the terms 


I<-> 


as arising from the inhomogeneity of turbulence. 

Unlike the non— scaled form, the energy equations of the macro-group contain 
the mode-couplings by the presence of the transfer functions. We also see the 
merit in these equations of being able to determine the spectral densities 
f and from the scaled energies 

0 

which are in reality the cumulative spectral distributions. 


C. Boundary Layer 

In boundary layers, the parameters (15) are restricted to the components 


l\,j 'JaT , with ^ ^ J 


( 17 ) 
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so that the energy equations (12) and (13) become reduced into the following 
particular form: 



~^o ,0 

K ^ - 

— ■ 1 O 
' ! 1 

(18a) 


- 

tz 

^ 1 • ^ 
'll 

(18b) 

1 1 (y-'i ') ^ 


V, 

(18c) 


X, Z 


(18d) 


The dissipation functions are negligible for large scale turbulence. 


D. Kinetic Representation 

It can be shown that the hydrodynamical equations of turbulence (la) and 
(lb) are simply the zeroth and first moments of the master equation 


(it V-F -V 





with the equivalence relation 



f V ~ 


(19) 


( 20 ) 


in the micro-dynamical state. The master equation has the advantage of being 
homogeneous, and its advective velocity does not cause a nonlinearity since V Is 
an independent variable. 

Like the hydrodynamical equations, the master equation can be submitted 
to group— scalings, so that the evolution of ^ ^ shown to be 

controlled by the transport property. 

£> y\(f z'ftjy . (21) 
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called diffusivity, with the evolution operator . This transport 

property attains its equilibrium hy a relaxation process, that is prescribed 
by A(t,t--o) and achieved by the diffusion of the trajectory at a diffusivity 
■ I f the two diffusivities and are structurally similar, a 

closure ±§ found. 

The analysis of the transport functions (18) will encounter many difficulties 
from the hydrodynamical approach. It is hoped that the above relaxational 
closure and the advantages of dealing with a homogeneous master equation of 
lesser nonlinearity in the kinetic approach will provide with a more powerful 
tool €or the determination of the transport functions. 


III. KINETIC RXJNDAim OF THE PRESSURE-^STRAIN CORRELATION 

Upon scaling the master equation to the macro-group, we have 

with the solution 

^ -A "Jclx A [t i-'c) ^ f _ 

''o 

By taking the moment and differentiate, we find 



( 22 ) 


(23) 


(24) 


(25a) 


i 
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(25b) 


and the divergence-free condition: 


v-u = 0. 


The'differentiation by the operator ^ has yielded the two contributions (25a) 
and (25b), and we have approximated the exact evolution operator /\. by the 
approximate operator XL/ , called the average propagator. Ncxv we rewrite 
these two terms in Fourier space, as follows: 


f 

p 1/ 0 


(27a) 


(2Zb) 




(28a) 


is an orbit function due to the unperturbed streaming, as opposed to the orbit 


function • \ f/'i r . 

. (j- **'> 

due to the perturbed path %(x^ by turbulence. 


(28b) 


The V-integration by parts gives 


(29a) 


(29b) 


Note that the effect of the orbit function! is small as compared to that of 

V 

i . The factor f(iyh as the role of securing the divergence- free condition. 


J 
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After the v- integration, we reduce (27) into: 



By definition (3a), the pressure 


Ch ^ 


\>Cl) =j4h pLkJ 

has a Fourier component 

k'- ^ - 

e/v] 


(30a) 


(30b) 


(31) 


so that the pressure- strain correlation becomes 



J 
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tensors from the cross- correlations. The. traces of the tensors do not possess 
indices. The diffusivities can play the role of integral operators in the 
integrations with respect to ~0 and . Note that (29) and (33) satisfy 
the divergence- free condition (2&), and that 

(35) 


is a boundary layer turbulence with the parameters (17). It leaves us the 


corre 


lation (33b). Here the reduction of the diffusivity 


K. 


by (34a ) has ,the benefit of confining ourselves to the velocity f luctuationf? 

, by avoiding the pressure- field fluctuations. 

" o 

The diffusivity from the auto- correlation has been investigated earlier by 
2 

Tchen . It is found that for a strongly stable boundary layer, i.e., 


va I /Y (36) 

the effect from ^ is negligible as compared to that from ^ . Here N is the 
Brunt-Vaisala frequency. The latter diffusivity can be investigated by the 
same method indicated above. Without going into the details of the calculation 
which will be reported at a later opportunity, we expect that 
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^ (V \vui K 

Z Z, MW 


( 37 ) 


The cascade transfer function is 


from our earlier work. The eddy viscosities are and f\ * 
spectral balance (18b) becomes 

I ■ 


(38) 


Hence the 


(39) 


For a steady state we differentiate (39) with respect to k, and obtain 


or 


by writing 


h/l 7 ul K°~ 

I 

((7uv>=£rW k'’F(k'j 

vn 


The spectrum is thus found to be 


where 


iJi) = ‘T Njvu j k 
< -ikKVK' 


~3 


(40a) 


(40b) 


(41) 


(42) 


(43) 


is a factor without dimension. The existence of the spectrum of power law 
can be seen from the experiments in Fig. 1. 'flrln spc'ctrum differs from the 
power law 

/i-\ 

(44) 




suggested in the literature. The two spectra (42) and (44) are in the ratio 


fe I /hi < I 


(45) 


i 
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stable turbulence. 



5-14 


for stable boundary layer. This small ratio is also verified by experiments. 

For a very stable boundary layer, where the conditon (45) is satisfied, the 
spectral distribution r;4) is similarly dominated by the balance between 
the redistribution and the cascade transfer, so that the power law is again 

valid. With the strong buoyancy, the spectrum /^j^^decays rapidly by going to 
large ^ , without giving the opportunity for efficiently. 
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SECTION 6 


Kinetic Equation of Turbulence 


ABSTRACT 

The different types of turbulence in applications to the atmospheric 
and oceanic motions, the propagation of light, the solitary waves, and the 
plasmas have different governing equations. W shall bring them to a 
canrroi Liouville form. In this manner, we hope that a single statistical 
method can be made available to treat the different types of turbulence. 
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I. FLUID EQUATIONS OF TURBULENCE 


Many basic phenomena in space, astrophysics, optics, atmospheric and 
oceanic sciences are nonlinear and random, i. e, they are in a turbulent 
state. A proper understanding of these phenomena depends critically on 
our ability to analyze the turbulent behavior. Although turbulence 
appear in different forms and with different basic equations for the 
description of their microdynamical-states, there exists certain similarity. 
To exploit the analogy in the attempt of finding a conmai form of equation 
of microdynamical state, we first consider incompressible and homogeneous 
turbulence, and write the governing Navier-Stokes equations, as follows: 




VU, =■ 0 


(ID 


( 2 ) 


Here is the velocity, P is the eonstant density 

the incompressible fluid, is the pressure, whose gradient enters into 

I ^ 

the field t. The variable in the symbol ( ) represents the total motion 

which can be decomposed into a mean and a fluctuation. The kinematic viscos- 
ity v’ , which is a transport coefficient from the molecular dissipation, is 
negligibly small compared with its turbulent counter part, and can be 
omitted in the analysis of the turbulent transport processes. The force of 
gravity forms a buoyancy, and together with the Coriolis term may be included 
in the ^ field. 

applying the divergence-free condition ( 2 ) to the equation of momentum 

n 

(1), ve obtain the following relation between C and 6i , in the form 

'V vwi 

V- E S ^ C3) 

vV ••4 tv w ' 
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called the "equation of state". 

The equations of motion of an inviscid ccirpressible fluid are as 


follows: 




f W. = f (4) 

■ ^ ^ t! ^ ^ 

-f f u. = 0 . (5) 

^ differentiating with respect to t and eliminating , 

we attain the relation 

^ C. . ^ ^ /I \ 

■I- (« 

The diffusion of a scalar A ] with chemical reaction E" is governed 

W 

by the following equation 

=E. v-^=o. (7, 

vAiere y? is the diffusion coefficient, and LC is the velocity field. 

Th® aaaa equations (7) are valid for the two-dimensional geostraphic 

A 

tuzioulence driven by a random field E . The geostraphic turbulence has 

-V* 

a voudzicity ^ 

^ 7X U , 

with Z'(Yh'X^) in two-dimensions. 

The equation of prapagatian of laser light has the form: 

Zik~ f -f- k^(t-i - = O . C83 

}' 

Here li, is the light field propagation in the z-diu:ection with an optical 
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wavenumber^ in a medium of fluctuating dielectric coefficient 2 . The 

term - represents the effect of the strong light intensity on the 


dielectric medium, causing the self-focusing of light. 

The equation (8) is called the "nonlinear Schrddinger equation 


1,2 
M ' 


It also applies to the solitons in plasmas, the atmosphere and the oceans. 
Here jg the envelop of the waves. 


The linearized form 


i- Vj^Ct+ 

^ 3 4 

IS called the parabolic equation of propagation. ' 




7>x^ 


( 9 ) 


Although the equation (9) is linear in fb , the product of two random 

A /\ 

functions lu makes the problem statistically nonlinear. In the form (9), 
the propagation of light resembles the diffusion (7) of a scalar. 


Recently a great deal of effort on the nonlinear wave phenomena has 

5 

been made by using the Korteweg-de Vries (KdV) equation: 

(\ -h LJ U. = 0, uHt/, L s u.— f A— . (10) 

where /A is the coefficent of dispersion. It has been shown that this 
adequately describes the long-time asymptotic behavior of wave motions in 
nonlinear dispersive systems, and that this asymptotic form can be 
derived either from the nonlinear Schrodinger equation (8) or from the 
perturbations of the equations of wave motions. 


1 
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II. MASTER EQUATION FOR THE DESCRIPTION OF THE MICRODYNAMICAL STATE OFTORBULENCE 


Although the fluid turbulence obeys different forms of differential 
equations (1) " (7) for the description of the microdynamical state, their 
kinetic representation can be cast into a unified master squation, or 
Klimontovich equation: 


+ L.) ~ ^ * 

Here ) is the distribution of the random variable ^ at time t~ 

and 'position X , and is a differential operator which takes 

HI tv ««« 

varied forms depending on the types of turbulence in consideration . 
These forms are found in the following. 

(i) For the inviscid incompressible turbulence as governed by (1) 
(3), ve have the differential operator 

^ H. sy.y-tE'^ / 


( 11 ) 


( 12 ) 


the equation of state 


A 

E 




(13) 


and the source 


- W K F' . 

^ _ ««v vn. 


(14) 


The limits of integration are understood to extend to the whole 
available domain, unless otherwise specified. The equivalence between the 
two representations, i.e. the fluid representation (1) - (3) on the one 
hand and the kinetic representation (11) on the other, can be obtained by 
writing 


f[v - U^(t,x)J • 


i 


(15) 
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With this equivalence, it is not difficult to transform the equation of 
state from the kinetic form (13) into the hydrodynamic form (3), and the 
master equation (11) into the hydrodynamic system (1) and (2) hy means 
of moments. 

(ii) For the compressihle turbulence, as governed hy (4) " (6), we 
have the same differential operator (12) and the same equation of state (13). 


The equivalence between the two representations is obtained by writing 
the distribution as 


II 

> > 

^ U(t,x)] , 

**• Iv 

(16) 

the equation of state as 



IfB = 


(17a) 

and the source as 

1/ j = 

— V-f k' 

(17b) 


By the equivalence relation (16), it is not difficult to transform 
the kinetic representation with formulas (11) and (16) into the fluid 
representation with formulas (4) " (6). 


(iii) For the diffusion and the two-dimensional geostrophic turbulence, 
as governed by (7), the Liouville equation (11) remains valid, with a 
differential operator 

I = i-v - y r' , 

and an equivalence condition 

v) = . (19) 

I *»«• 

(iv) For the light propagation, we can introduce 


3 
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and transform the Schrodinger equations (8) and (9) into 

L ^ ^j] • 

with 

f ^-ivk"(i-i-iii^r) ( 22 .) 

in nonlinear propagation, and 

lj> = - (s - I ) (22b) 

in 1 in ear propagation. 

Note that (21) is already in a form analogous to the master equation 
so that further transformation into the kinetic representation becomes 
unnecessary. The same argument holds for the KdV equation (10). 


(v) The plasma turbulence is governed by (11), now called the Vlasov 
equation, with the differential operator (12). The E - field represents 
the electrostatic field with unit charge and mass, The equation of state in 
the form (13) remains valid, with the source 


•V* W' 


(23) 


where is the constant plasma frequency. The distribution function has 
the normalization condition 




f = 

where *^he number-density. 


(24) 
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TYPES OF TURBULENCE LIOUVILLE DISTRIBUTION DIFFERENTIAL SOURCE 

EQUATION OPERATOR 
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In Table 1, we have listed the different types of turbulence. Their 
governing equations are written in the master form with appropriate 
differential operators. The self-consistent field is related to the 
distribution function by the equation of state in the form (13) through 
a source A . This reduction of several types of turbulence into a unified 
master equation suggests that they can all be treated by a single 
statistical method. 


The description of the microdynamical state of turbulence by an equation 
of the master form presents the advantages that the equation is homogeneous, 
and that the replacement of the velocity function by the independent 

variable V prevents any terms connected kith the velocity from becoming 
nonlinear, as found in the differential, operators L and the sources (a) " (c) . 


The only surviving nonlinear term arises from the E -field. 

Fortunately this does not require our immediate attention, since it can be 


provisionally considered as a given random force until after the closure. 


A 

Thereafter, E is determined by the equation of state with the proper source. 

This nonlinear term describes the wave-particle interation 
explicitly, and therefore will enable the derivation of the nonlinear 
Landau damping in fluids in a direct manner. The fluid representation may 
describe such an interaction only indirectly through a series of moments. 
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3. Kinetic Hierachy of Turbulence 


In the following lines we shall be concerned with the incompressible 
Navier-Stakes turbulence. The microdynamical state of turbulence is 
described in the kinetic representation by the master equation 



(25) 


with a differential operator 

t = VF - V7* V- E.‘l 

and a self-consistent field 

! -rJ ' 

so that the pressure |? satisfies the Poisson equation 

vi-rz-i 


(26) 


(27) 


(28) 


by (13). The solution can be written in the form: 




(29) 


W calculate 


--V jk/dv '. — L 

•“ ‘V - - fTT/x-x'l / — 


(30) 


by differentiating (29), and obtain 

A 




(31) 
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by introducing the integral operator 



( 32 ) 


The term (31)yields a second order nonlinearity in the master equation. 


The master equation, which describes the microdynamical state of 
turbulence, contains fluctuations with all minute details which are 
unnecessary, if not impossible, in a statistical study. A coarse-graining 
procedure is the "global ensemble average" with the operator 

^r<r V . 

Such an average is deterministic and may vary with or ^ ^ . The 

deviation from this average is obtained by the operator of fluctuation 

-A/ " — 

/\ =• t-A , 

where "1" is the unit operator . 

applying the operators ^ ,A to (25), we obtain: 

+Z)J = - A ^ C , ( 33 ) 


and 


or 




(34) 


(35) 


The operators operate on the functions which follow. 

The equation (33) describes the evolution of the distribution 
in a turbulent meciium which presents a turbulent collision 
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representing the statistical effects of the fluctuations upon 
that the collision, as defined hy 


Note 



(36) 


is the derivative of the Reynolds stress 



(37) 


in the phase space. In this sense the equation (33) can he called the 
Reynolds equation in the phase space, and the equation of. evolution 
of the fluctuation, (34) or (35), is the Friedmann equation, as obtained 
hy the Reynolds decomposition 


I = i-A • 


(38) 


1^ this decomposition, we can separate (32) into two components, as 
follows : 

“I'M “ ' f y 

i-}f = 


(39b) 


where 


= d f{t xj A 


(40) 


is the p^ir distribution function. Upon substituting (39) into the 
ensemble average of (32) and subsequently into (33) , we get the following 
equation of evolution of ^ : 


J 



6-13 


! (y, V-/ <; v) j f ft*,v '^1 + C , (41) 


with 




(42) 


The dependance of the collision upon the pair distribution function 

is expected. If we consider a source of the type (23), ^\e reduce (41) to 

the BBGKY (Bogoliubov, Bom, Green, Kirkwood and Yvon) hierarchy of plasmas . 

7-12 

Our general form (41) applies to other types of turbulence too. 

The equation of evolution of ^ in (41) will be called the kinetic 
equation of turbulence. remains to determine the collision as an ex plicit 
function of , i.e. in the form: 


C = S(/| 


(43) 


where is a collision operator. 


1 
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